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1  INTRODUCTION 


The  modeling  of  wind  flow  fields  over  complex  terrain  is  a  difficult  problem 
that  is  not  new  to  the  field  of  Meteorology  (Queney,  1948;  Liu  and  Goodin,  1976). 
The  complicated  flow  patterns  produced  from  air  flowing  past  irregular  surfaces 
and  obstacles  are  similar  to  wave  and  eddy  features  seen  in  rushing  stream  wa¬ 
ter.  Calculations  necessary  to  predict  such  flow  patterns  are,  needless  to  say, 
computationaly  intensive  (Deardorff  et  al.,  1984;  Endlich,  1984).  To  ease  the 
computational  burden,  model  developers  make  many  simplifying  assumptions 
regarding  flow  characteristics  and  boundary  conditions  (Kao,  1981).  Unfortu¬ 
nately  this  often  leads  to  an  over  simplified  flow  solution.  For  example,  even  low 
hills  can  significantly  affect  wind  flow  patterns  (Smedman  and  Bergstrom,  1984). 
Nastrom  et  al.  (1987)  note  that  complex  terrain  induces  significant  variability 
to  local  wind  fields. 

Real  atmospheric  flow  patterns  are  marked  by  complicated  trajectories  in  3- 
dimensional  space  (Endlich,  1967).  Therefore,  if  realistic  wind  fields  over  complex 
terrain  are  to  be  realized  in  model  outputs,  solutions  must  integrate  horizontal 
flow  behaviors  v  i;h  those  in  the  vertical.  When  the  3-dimensional  mechanical 
effects  of  terrain  on  flow  fields  are  not  considered,  then  high  speed  flows  through 
mountain  passes  and  deviations  in  flow  direction  along  the  edges  of  solid  obsta¬ 
cles  will  not  be  correctly  modeled  (Sherman,  1978). 

1.1  Literature  Review 

Why  is  so  much  effort  expended  in  developing  models  to  simulate  these  compli¬ 
cated  wind  field  behaviors?  Accurate  representation  of  wind  flow  fields  is  im¬ 
portant  to  realistically  determine  where,  and  in  what  quantities,  air  pollutants 
may  be  found  (Davis  et  al.,  1984;  King  and  Bunker,  1984).  This  also  applies 
to  the  simulation  of  dust  and  smoke  aerosol  models  used  to  simulate  battlefield 
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obscurants  (Hoook  et  aJ.,  1987).  The  prediction  of  low  level  aviation  hazard  ar¬ 
eas  in  complex  terrain  (Lee  and  Hansen,  1989),  or  deciding  where  to  place  wind 
power  stations  (Bhumralkar  et  al.,  1980)  are  other  applications  areas  that  require 
detailed  wind  field  predictions.  Because  it  is  not  practical  to  measure  wind  vec¬ 
tor  field  properties  at  all  locations  over  a  selected  area,  wind  modeis  have  been 
developed  to  estimate  winds  at  other  locations  using  only  a  selected  amount  of 
available  data  (Bhumralkar  et  al.,  1980). 

Wind  flow  properties  may  be  studied  at  many  different  levels  of  detail,  and 
that  presents  design  problems  to  model  developers.  A  major  issue  among  atmo¬ 
spheric  modelers  has  been  agreeing  upon  the  scale  employed  to  resolve  flow  field 
features  (Asculai  et  al.,  1984;  Bannon  and  Zehnder,  1989).  Naturally,  many  mod¬ 
elers  would  prefer  to  study  flows  in  great  detail,  but  this  has  not  always  been 
practical  because  of  the  computational  burdens  that  detailed  models  present. 
However,  each  major  improvement  in  computer  storage  capacity  and  speed  eases 
this  problem  to  some  degree.  The  trend  is,  therefore,  to  model  flow  fields  us¬ 
ing  ever  smaller  grid  spacing  -  from  100’s  of  kilometers  to  less  than  1  kilometer 
(Erasmus,  1986;  Fast  and  Takle,  1988).  Model  performance  in  this  proposal  will 
be  examined  in  a  high  resolution  grid  (e.g.,  <  1  km  grid  spacing). 

Consistent  approaches  that  have  become  widely  accepted  in  the  field  of  com¬ 
plex  terrain  wind  flow  modeling  include  the  use  of  sigma,  or  terrain  following 
coordinate  systems  (Mahrer  and  Pielke,  1975;  Mass  and  Dempsey,  1985),  and 
the  assumption  of  neutral  stability  throughout  the  modeled  flow  field  (Beljaars 
et  al.,  1987).  Also,  today’s  modelers  seem  to  be  placing  more  emphasis  on  the 
application  of  dynamic  flow  theory  (Beljaars  et  al.,  1987;  Smith,  1988;  Bannon 
and  Zehnder,  1989)  than  on  empirical  parameterization  of  flows  (Han  et  al.,  1982; 
Marwitz,  1983;  Dinar,  1984;  Lee  and  Kau,  1984).  All  these  approaches  were  in¬ 
tegrated  to  a  significant  degree  in  developing  the  model  to  be  evaluated  in  this 
proposal  (Ludwig  and  Endlich,  1988). 
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In  the  relatively  short  history  of  wind  modeling,  some  model  developers  took 
more  complex  3-dimensional  approaches  than  others,  to  develop  very  accurate 
representations  of  real  flow  field  properties  (Dickerson,  1978;  Sherman,  1978). 
For  example,  Ludwig  et  al.  (1990)  developed  a  model  that  budgets  gravita¬ 
tional  potential  energy,  mgy,  and  kinetic  flow  energy,  jmt/a,  parallel  to  wind 
field  streamlines  where  m  is  mass,  g  is  acceleration  due  to  gravity,  y  is  distance 
along  the  vertical  axis,  and  v  is  the  flow  velocity  (Sears  et  al.,  1982).  This  bud¬ 
geting  process  allows  for  lateral  mass  transfer  determinations  to  be  made  from 
mechanical,  terrain  forcing  effects.  Lateral  mass  transfer  processes  are  compared 
to  the  overly  simplified,  vertical  forcing  method  as  shown  in  Figure  1.  The  part 
of  the  lateral  transfer  potential  model  that  actually  calculates  wind  fields  is  orig¬ 
inally  from  a  previous  endeavor  (Ludwig  and  Endlich,  1988)  and  is  referred  to 
as  the  Winds  on  Critical  Streamline  Surfaces  (WOCSS)  code.  There  was  also  a 
precursor  to  the  WOCSS  model,  described  by  Ludwig  et  al.  (1986),  which  incor¬ 
porated  an  initialization  scheme  that  used  a  very  small  amount  of  observation 
data. 

Figure  2  illustrates  that  the  WOCSS  code  evolved  from  well  researched  and 
widely  accepted  wind  modeling  principles.  Consideration  of  lateral  mass  transfer, 
although  a  unique  aspect  of  the  WOCSS  code,  is  not  the  only  strong  point  of  this 
particular  model;  the  WOCSS  model  is  diagnostic,  which  means  that  it  does  not 
require  time-dependent  boundary  conditions  (Endlich,  1984).  Prognostic  models 
are  used  to  make  time-dependent  forecasts  via  the  integration  of  conservation 
relations,  which  involves  many  added  calculations;  diagnostic  models,  however, 
“are  very  economical  and  appear  to  be  effective  mesoscale  tools  when  the  dom¬ 
inant  forcing  is  terrain  according  to  Pielke  (1984,  p.480).  As  a  result  of 
selecting  the  diagnostic  approach,  Ludwig  et  ai.  (1990)  were  able  to  develop 
a  compact  code  that  could  be  installed  on  micro-computers,  and  could  also  be 
executed  relatively  faster  than  many  other  models  that  simulate  wind  flow  over 
rough  terrain. 
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LATERAL  FORCING  CASE 


Figure  1.  Vertical  forcing  wind  flow  model  solutions  (a)  do  not  make  the 
necessary  potential  energy  considerations  to  simulate  the  effect 
of  lateral  mass  transfer  that  occurs  in  nature  (b). 


Figure  2.  This  figure  illustrates  a  map  of  the  reference  links,  or  information  scheme,  of 
central  importance  to  this  research  work.  Ludwig  and  Endlich  (1988)  describe  the  wind 
model  selected  for  validation  In  this  study.  The  foundation  for  the  design  of  the  model  is 
based  on  references  to  well  known  work  in  the  field  of  wind  modeling.  Foundations  to 
some  of  Henmi’s  (1988)  modeling  work,  parallel  to  Ludwig’s  work,  are  also  illustrated  in 
this  figure. 
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1.2  Hypotheses 

The  WOCSS  complex  terrain  wind  flow  model  is  designed  to  efficiently  simulate: 
1)  lateral  mechanical  speed  up  effects  of  air  flows;  and  2)  changing  wind  directions 
along  elevated,  solid  boundaries.  These  considerations  should  make  the  WOCSS 
model  outputs  significantly  more  reliable  than  other  models  that  do  not  address 
lateral  forcing  in  complex  terrain.  Therefore,  the  major  assumption  made  prior 
to  the  commencement  of  this  research  work  was  that  the  WOCSS  model  was 
significantly  more  reliable  than  basic,  vertical  forcing  wind  models  because  the 
WOCSS  model  makes  considerations  for  lateral  flow  adjustment;  it  was  further 
assumed  that  this  increase  in  reliability  was  not  defeated  by  excessive  additional 
processing  time. 

1.3  Research  Objectives 

The  WOCSS  model  has  not  yet  been  independently  evaluated,  and  Ludwig  et 
al.  (1990)  admit  that  more  testing  and  evaluation  of  the  WOCSS  model  are 
required  to  validate  the  hypothetical  strengths  of  their  model.  Therefore,  the 
primary  objective  of  this  research  work  was  to  implement  the  WOCSS  model 
on  selected  computer  resources,  and  establish  quantitative,  statistical  measures 
of  the  WOCSS  complex  terrain  model’s  performance  characteristics  via  compar¬ 
isons  made  using  high  resolution  meteorological  wind  data  sets  (both  surface 
and  upper-air  wind  profiles),  with  a  basic  vertical  forcing  wind  model  developed 
by  Henmi  and  Tabor  (1988).  Henmi’s  model,  used  in  this  research  work,  em¬ 
ployed  an  objective  wind  analysis  approach  and  will  hereafter  be  referred  to  as 
the  OBWIND  model. 

The  WOCSS  evaluation  process  was  carried  out  using  three  separate  com¬ 
plex  terrain  surfaces.  Finally,  the  WOCSS  model  was  evaluated  using  a  remote, 
fourth  complex  terrain  surface  under  sparse  meteorological  data  initialization 
conditions;  this  study  was  carried  out  to  further  appraise  the  potential  utility 
of  the  WOCSS  model  for  use  in  supporting  U.S.  Army  research  of  battalion- 
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regiment  size  battlefield  problems.  Figure  3  is  a  conceptual  illustration  of  the 
basic  study  design. 

The  WOCSS  model  evaluation  process  carried  out  in  this  study  represents 
necessary  third  party  analysis  as  a  step  toward  model  validation.  Simply  defin¬ 
ing  model  validation  has  become  a  significant  exercise  in  itself.  However,  it  is 
generally  considered  by  many  scientific  communities  to  be  a  process  leading  to 
the  determination  of  whether  a  model  reliably  predicts  events  that  actually  occur 
in  the  real  world  (Hodges  and  Dewar,  1990).  The  comparisons  described  in  the 
previous  paragraph  were  developed  to  identify  measures  of  the  WOCSS  model’s 
ability  to  predict  atmospheric  flow  structures,  as  they  actually  occur  in  nature; 
and  the  contrast  of  model  types  used  in  making  these  comparisons  was  also  in¬ 
tended  to  improve  physical  insight  into  fundamental  properties  of  atmospheric 
motions  in  complex  terrain  (Holton,  1979). 


Figure  3.  Basic  wind  model  evaluation  study  design. 
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2  METHODS 


The  steps  taken  to  develop  this  study  included  the  selection  of  terrain  and 
meteorological  data,  and  preparation  of  the  WOCSS  and  OBWIND  models  for 
case  study  runs. 

2.1  Terrain  Data  Bases 

Four  complex  terrain  sites  were  selected  for  use  in  this  validation  study.  All 
four  complex  terrain  areas  were  set  to  20.0  km  x  20.0  km  sizes,  representative 
of  selected  battalion-regiment  size  battlefield  engagement  areas.  As  a  result, 
these  four  areas  were  immediately  representative  of  mesoscale1  dimensions,  in 
meteorological  terms. 

The  four  sites  selected  were  1)  the  Orogrande  region,  located  on  the  White 
Sands  Missile  Range  (WSMR)  in  the  Tularosa  Basin,  south  central  New  Mexico; 
2)  the  National  Training  Center  (NTC),  Fort  Irwin  California,  which  is  located 
south  of  Death  Valley;  3)  the  Sierra  Nevada  foothills,  near  Red  Bluff,  Califor¬ 
nia;  and  4)  a  complex  terrain  region  approximately  25  km  southwest  of  Pilsen, 
Czechoslovakia.  The  general  locations,  and  a  contoured,  topographic  plot,  of 
sites  1-4  are  listed  in  Table  1  and  illustrated  in  Figs.  4-7. 

Mesoscale  meteorological  observation  programs  have  been  conducted  on  com¬ 
plex  terrain  sites  1-3.  In  fact,  one  of  the  most  extensive  and  detailed  meteorolog¬ 
ical  observation  programs  that  has  ever  been  conducted,  took  place  in  the  site  3 
region;  this  program  was  called  Project  WIND  (Cionco,  1989).  Therefore,  site 
3  will  be  referred  to  hereafter  in  this  report  as  the  Project  Wind  site.  Project 
WIND  will  be  referred  to  later  when  meteorological  data  is  discussed. 

‘Pielke  (1984,  p.l)  defines  mesoscale  as  a  horisontal  scale  “on  the  order  of  a  few  kilometers 
to  several  hundred  kilometers  or  so,  with  a  time  scale  of  about  1  to  12h." 
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General  location  of  the  NTC  site,  marked  by  the  small  box  in  southeastern  California  (a);  and  a  topographic 
a  20  X  20  km  site  area  (b),  where  each  grid  cell  represents  a  1  X 1  km  area,  with  a  contour  interval  of  50  m. 
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Figure  6.  General  location  of  the  Project  WIND  site,  marked  by  the  small  box  in  northeastern  California  (a);  and  a  topographic 
map  of  the  20  X  20  km  site  area  (b),  where  each  grid  ceil  represents  a  1  X  1  km  area,  with  a  contour  interval  of  25  m. 
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Figure  7.  General  location  of  the  Czechoslovakian  site,  marked  by  the  small  box  near  the  West  German  border  (a); 
and  a  topographic  map  of  the  20  X  20  km  site  area  (b),  each  grid  cell  represents  1  X  1  km,  with  a  25  m  contour  interval. 


Site  4  was  selected,  in  contrast  to  the  Project  WIND  site,  as  a  remote,  com¬ 
plex  terrain  area  to  test  under  sparse  data  conditions;  no  actual  meteorological 
observation  data  was  available  for  site  4  in  this  study.  This  exemplifies  the  cir¬ 
cumstances  under  which  the  WOCSS  model  could  be  used  for  support  purposes 
under  extenuating  field  conditions. 


Table  1.  Approximate  Latitude/Longitude  of  the  southwest  corner  of 
the  four  selected  20  km  x  20  km  terrain  sites. 


Site  No. 

Site  Name 

Latitude 

Longitude  | 

IHUH 

WSMR,  NM 

32°  1!)'  N 

1 06°  13'  W  | 

11 — 

NTC,  CA 

35°23'  N 

Project  WIND,  C A 

40°8'  N 

122° 18'  W 

iikhi 

Pilsen,  Czech. 

49°30'  N 

12° 10'  E 

Digital  terrain  elevation  -  level  1  data  (DTED)  was  acquired  for  sites  1-1  from 
the  Defense  Mapping  Agency  (DMA).  The  elevation  data  were  relative  to  mean 
sea  level  (MSL),  and  for  the  purposes  of  this  research  work,  only  elevation  data 
records  were  employed.  DTED  -  level  1  elevation  values  are  given  in  meters.  The 
DTED  data  were  provided  in  1°  x  lc  files  that  were  identified  bv  their  southwest 
corner  coordinates.  The  lat.it  ude/longit  ude,  or  A .r.  Ay,  grid  spacing  of  DTED  - 
level  1  elevation  data  is  3  x  3  arc  seconds.  Since  there  are  3600"  in  1°  of  longitude 
or  latitude  and  111000  m  ~  1°  of  latitude  or  longitude,  both  the  Ar.  Ay  interval 
is  yE—flllOOO)  ni  ~  92.50  in,  which  led  to  the  establishment  of  100  m  as  the 
highest,  resolution  possible  in  the  x,y  terrain  grid  spaces. 


13 


One  hundred-meter  grid  spacing  is  an  unusually  dense  meteorological  grid 
and  will  provide  a  degree  of  resolution  in  flow  field  simulations  that  has  been 
difficult  to  attain  in  previous  modeling  work,  basically  due  to  hardware  related 
computational  limitations  and  meager  availability  of  mesoscale  meteorological 
data  bases.  Orlanski  (1975)  discussed  several  resolution  planes,  including  100 
m,  that  can  be  used  to  identify  gravity  wave  features2  that  can  contribute  sig¬ 
nificantly  to  ambient  atmospheric  properties,  such  as  the  degree  of  turbulence 
present  at  a  particular  location.  Although  the  research  work  presented  here  is 
focused  more  on  lateral,  rather  than  on  vertical  gravity  wave  effects,  the  point  is 
that  100-m  grid  spacing  is  expected  to  resolve  meteorologically  significant  devi¬ 
ations  (wavelengths  =  A)  in  simulated  flow  fields. 

Hansen  and  Flanigan  (1989)  developed  a  program  to  generate  modified  terrain 
elevation  data  in  universal  tranverse  mercator  coordinates  from  one  to  four  1° 
DTED-level  1  data  cells,  and  this  program  was  used  to  develop  the  100-m  grid 
systems,  as  well  as  two  additional,  gradually  coarser,  Ax,  A y  grids  for  use  in  this 
validation  study  (the  selection  was  restricted  to  this  number  of  cases  because  of 
time  limitations). 

The  purpose  in  selecting  more  than  one  Ax,  Ay  grid  mesh  was  to  examine  the 
effect  of  varying  grid  spacing  on  the  wind  models’  abilities  to  resolve  terrain  effects 
in  their  atmospheric  flow  structure  solutions.  Since  the  highest  resolution  limit 
was  set  by  the  DTED  spacing  format,  the  problem  was  to  establish  a  reasonable 
low  resolution  Ax,  Ay  limit,  as  well  as  an  intermediate  Ax,  Ay  resolution.  Several 
factors  led  to  the  selection  of  the  low  Ax,  Ay  grid  resolution  limit. 

2Hooke(1986)  describes  gravity  waves  as  mesoscale,  incompressible  environment  responses 
in  the  vertical.  These  features  are  coupled  to  the  effects  of  local  gravitational  force  and  the 
degree  of  stable  stratification  in  the  local  atmospheric  column. 
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One  of  the  primary  factors  in  this  selection  was  directly  related  to  the  20.0 
km  x  20.0  km  size  of  the  terrain  surfaces  themselves.  This  lower  limit,  or  coarse, 
resolution  grid  system  should  contain  enough  z,  y  data  points  to  depict  reasonable 
spatial  averages  of  the  associated  wind  flow  field  properties  (Brown,  1991). 

Smaller  scale  perturbations,  induced  by  complex  terrain  in  the  basic  flow 
field,  u,  would  be  of  more  interest  at  higher  Az,  Ay  resolutions  (Holton,  1979), 
such  as  at  the  Az  =  100-m  resolution.  Other  factors  leading  to  the  selection  of 
the  coarse  resolution  grid  spacing  included  a  much  more  detailed  consideration 
of  perturbation  properties  in  simulated  flow  structure  solutions. 

A  minimum  of  3  grid  points  are  required  to  fully  describe  a  complete  sinusoidal 
wavelength  feature  (or  perturbation);  this  is  equivalent  to  a  distance  of  2Az,  or 
2Ay  (see  Fig.  9),  and  is  referred  to  as  the  resolution  limit,  /Z,  of  an  z,y  grid 
(Riegel,  1974).  R  can  be  expressed  as 

(1) 

where  pi  =  Here  A  is  the  wavelength  of  the  perturbation  of  interest,  and  Az 
is  the  grid  spacing  in  the  z  coordinate  direction  (Riegel,  1974). 


Wavelength  of 
Perturbation  or 
Displacement 
Feature 


Figure  8.  Minimum  number  ol  grid  points  required  to  resolve  a  dsplacement  feature. 
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When  A  =  2Ax,  then  R  =  0,  perturbation  features  with  A  <  2Ax  will  not 
be  resolved  by  the  grid.  The  error,  in  percent,  that  can  be  made  by  incorrectly 
phasing  a  perturbation  <  A  through  a  grid  can  be  expressed  as  (1  -  R)  •  100. 

Examination  of  elevated  terrain  features,  on  1:250,000  scale  maps,  in  the 
selected  WSMR  and  NTC  terrain  areas  revealed  that  significant  terrain  effects 
on  atmospheric  flow  properties  would  probably  begin  to  disappear  at  A  >  2000 
m.  A  2000  m  A  represents  10%  of  the  x  or  y  linear  limit  (20  km).  Since  2  Ax  is 
♦he  perturbation  resolution  limit  of  a  grid  network,  for  a  perturbation  A  =  2000 
m,  Ax  must  <  1000  m.  Therefore,  1000  m  was  selected  as  the  lowest  resolution 
Ax,  Ay  grid  spacing. 

Table  2.  The  (1  —  /?)  •  100  resolution  limit  error,  in  percent,  that  is 
possible  in  phasing  perturbations  <  A  through  a  grid  with  a  given 
resolution  of  Ax. 


Ax  (meters) 

P 

3000 

erturb* 

2000 

jtiona 

1000 

<  A  (n 
500 

leters) 

250 

100 

^  3000  1 

100 

100 

100 

100 

100 

Em 

2000 

100 

100 

100 

100 

100 

100J 

1000 

58.7 

100 

100 

100 

100 

EEH 

250 

IKES 

10.0 

36.3 

100 

MM 

100 

100 

0.73 

1.6 

6.5 

mi 

76.6 

100 

Note  in  Table  2  that  for  Ax  =  1000  m,  the  resolution  error  (in  percent)  for 
all  listed  perturbations  with  a  A  <  2000  m  is  always  100%,  but  at  Ax  <  250  m 
there  are  significant  decreases  in  error  as  A  increases.  For  Ax  =  250  m,  error  is 
only  10  %  for  A  =  2000  m  and  is  futher  reduced  to  less  than  2%  when  Ax  =  100 
m.  There  are  also  other  reasons  to  support  the  Ax  =  250  m  spacing  other  than 
the  low  10%  error  threshold  it  has  for  A  <  2000  m. 
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Table  3  summarizes  the  basic  properties  of  three  20  x  20  km  grid  systems 
with  1)  1000;  2)  250;  and  3)  100  m  Ax,  Ay  grid  point  spacing,  respectively. 
Note  the  exponential  increase  that  occurs  in  the  52  of  all  points  column  when 
Ax,  Ay  spacing  is  changed  from  1000,  250,  and  100  m,  respectively.  This  growth 
behavior  is  illustrated  in  Fig.  9(a).  Figure  9(b)  is  a  log-linear  plot  of  the  same 
information,  but  it  is  closer  to  a  straight  line  trend  -  and,  therefore,  easier  to 
describe.  A  straight  linear  projection,  y  =  mz  +  b,  from  100  to  1000  m  produces 
a  mean,  intermediate  spacing  value  of  550  m  (see  Fig.  9(c)),  but  a  similar  log- 
linear  projection,  log  y  =  mx  +  b,  produces  a  mean,  intermediate  grid  spacing 
value  of  316.23  m  (see  Fig.  9(d)). 

Table  3.  The  basic  two  dimensional  grid  properties  of  a  20  x  20  km 
area  for  Ax  and  Ay  both,  respectively,  equal  to  1000,  250,  and  100 
m. 


Grid  System 

Ax,  Ay 

x  x  y  points  [| 

53  of  all  points 

1 

1,000  m 

21  x  21  1 

441 

2 

250  m 

81  x  81  | 

6,561 

3 

100  m 

201  x  201  || 

40,401 

The  nearest  multiple  of  the  upper  Ax  spacing  limit  (1000  m)  less  than  316.23 
m  is  250  m.  Reasons  to  select  a  multiple  of  1000  m  that  is  <  316.23  m  include: 
1)  that  it  facilitates  scaling  of  1000  -  100  m  Ax  cross  grid  system  results;  and  2) 
Table  2  suggests  the  use  of  a  smaller  Ax  to  reduce  error  in  A  resolution.  Because 
of  this,  and  the  exponential  nature  of  the  53  of  all  points  (see  Table  3),  Ax  =  250 
m  was  selected  as  the  intermediate  grid  resolution.  This  final  determination 
established  the  selected  set  of  Ax  resolutions  (1000,  250,  and  100  m)  that  were 
used  in  this  study.  Selection  of  this  set  of  Ax  spatial  resolutions  for  use  in  this 
study  required  further  precision  to  be  employed  in  defining  the  meteorological 
scales  of  motion  that  would  then  be  simulated  by  the  WOCSS  and  OBWIND 
models. 
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Linear  Scale  behavior  of  the  number  of  Log-Linear  behavior  of  the  number  of 

grid  points  selected  {2  grid  points  selected 
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Because  the  coarse,  Az  =  1000  m,  grid  was  established  to  resolve  A  >  2, 000 
m,  it  qualifies  as  a  mesoscale  mesh  (Pielke,  1984).  However,  Table  4  implies  that 
the  remaining  two  grid  structures  (Az  =  250,  and  100  m)  fall  into  a  different 
scale  -  the  microscale.  All  three  selected  Az  grid  systems  can  be  more  precisely 
defined  from  the  information  provided  in  Table  4.  For  example,  the  Az  =  1000 
m  grid  can  be  classified  as  a  meso  7  network,  because  its  A  resolution  limit  is 
2Az  =  2  km.  Similarly,  the  Az  =  250  and  100  m  grids  have  the  potential  to 
resolve  micro  a  scale  systems,  because  their  A  resolution  limits  are,  respectively, 
500  and  200  m. 


Table  4.  Selected  meteorological  spatial  (Az)  scale  definitions.  3 


Meteorological  Scales 

A 

Macro  a 

A  >  10,000  km  |; 

Macro  0 

<  A  < 

Meso  a 

HHSIBHi 

2,000  km  | 

Meso  0 

Meso  7 

2  km 

<  A  < 

Micro  a 

<  A  < 

—i 

Micro  0 

20  m 

<  A  < 

Micro  7 

A  <  20  m  | 

Table  5  lists  statistics  that  further  characterize  the  four  terrain  sites  selected 
for  use  in  this  study.  These  terrain  properties,  even  though  they  are  extremely 
significant,  are  not  really  results  themselves;  but  they  establish  the  context  within 
which  to  interpret  the  wind  model  test  results.  The  percentage  change  between 
250  -  100  m  Az  elevation  data  bases  were  not  significant,  but  some  1000  -  100  m 
changes  were  noteworthy;  therefore,  they  were  selected  for  inclusion  in  Table  5. 

3Orlanski  (1975). 
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The  percentage  change  between  all  the  1000  •  100  m  elevation  data  mean 
elevations  was  insignificant  (less  than  1%)  and  almost  normally  distributed  (see 
Table  5).  Minimum  elevations  and  distribution  frequencies,  inferred  from  the 
±1.0  standard  deviation  comparisons,  were  also  fairly  stable  between  1000  - 
100  m  (less  than  1%  and  3%,  respectively).  Both  the  mean  elevation  and  the 
±  1.0  standard  deviation,  at  100,  250,  and  1000  m  Ax  grid  spacing,  for  all  the 
terrain  data  bases  are  illustrated  in  Fig.  10.  Therefore,  variation  in  grid  spacing 
selection  between  100  -  1000  m  appears  to  have  little  relative  effect  in  changing 
each  site’s  mean  elevation  and  distribution  variance.  The  behavior  of  each  site’s 
mean  elevation  value  and  distribution  behavior,  therefore,  appears  stable,  in 
relative  terms,  when  the  grid  spacing  is  changed  (see  Fig.  10). 

Table  5.  Variations  in  the  terrain  elevation  data  bases  (m  above 
MSL),  as  a  function  of  grid  spacing,  Ax  (m),  for  the  listed  sites: 

1)  WSMR;  2)  NTC;  3)  Project  WIND;  and  4)  Czechoslovakia.  The 
variations  are  presented  in  terms  of  the  1000  m  - 100  m  Ax  differences. 

These  differences  are  then  expressed  as  percentage  change  between 
the  1000  m  -  100  m  Ax  terrain  data  bases. 


SITE 

Ax 

Terra 

MEAN 

in  Statistics  (inc 
MIN 

uding  lOOO-lOOm 
MAX 

A  %) 

STD.  DEV. 

1 

1000 

100 

1266.21 

1268.00  .1% 

1219.00 

1218.00  -.1% 

47.09 

45.81  -2.7% 

2 

EE3 

1036.84 

1043.81  .7% 

678.00 

677.00  -.2% 

■vTTv^pwn 

3 

m 

161.84 

160.68  -.7% 

83.00 

83.00  0% 

| 

69.56 

68.95  -0.9% 

■ 

1IIS 

490.48 

488.57  -.4% 

377.00 

374.00  -.8% 

| 

67.48 

65.76  -2.6% 

20 


o  o  o  o  o  o  o 

o  o  o  o  o  o 

CN  O  00  U3  CN 


|8A8|  688  U68UJ  8A0qB  UJ  U|  *Z  'U0flBA8|3 
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However,  comparisons  of  maximum  elevation  data  revealed  that  maximum 
elevations  could  change  as  much  as  +  31  %  (Table  5,  for  the  Project  WIND 
site)  when  Ax  grid  resolution  improved  by  a  factor  of  x  10  (e.g.,  the  1000  -  100 
m  comparison);  this  results  from  the  “smoothing"  effect  that  occurs  at  lower 
resolutions.  Figures  11-14  illustrate  this  “smoothing”  effect  when  Ax  resolution 
is  decreased  from  250  m  to  1000  m  for  the  WSMR,  NTC,  Project  WIND,  and 
Czechoslovakian  DTED  bases.  The  “smoothing”  process  detects  fewer  extreme 
elevation  points,  and  this  detection  failure  is  not  consistent  -  it  depends  on  how 
the  grid  mesh  is  registered  to  a  terrain  surface  and  will  also  vary  significantly 
between  different  terrain  surfaces.  From  Figs.  11-14  it  is  evident  that  careful  Ax 
selection,  for  the  purposes  of  studying  complex  terrain  effects,  is  critical.  Note 
that  the  NTC  site  has  the  greatest  variability  and  range  in  elevation,  while  the 
WSMR  site  elevation  has  the  smallest  extremes  and  variability.  Also,  the  distri¬ 
bution  of  elevation  data  for  the  Project  WIND  Site  produces  standard  deviations 
closely  resembling  those  of  the  Czechoslovakian  Site  (see  Table  5).  The  influence 
of  these  Ax  grid  variations  on  the  WOCSS  model  will  be  discussed  in  the  Results 
section. 


2.2  Meteorological  Data  Bases 

Project  WIND  was  a  cooperative  field  study  designed  to  collect  comprehensive 
meteorological  data  over  complex  terrain  for  meteorological  model  evaluation  and 
improvement  purposes  (Cionco,  1989).  Project  WIND  was  designed  to  address 
mesoscale  and  microscale  meteorological  domains  (see  Table  4).  The  meteorolog¬ 
ical  data  collected  during  Project  WIND  included  measurements  of  wind  speed 
and  direction,  turbulence,  temperature,  humidy,  pressure,  solar  radiation,  soil 
heat  flux,  precipitation,  upper  air  soundings,  and  satellite  images.  Observation 
sites  were  scattered  over  a  200  km  x  200  km  domain  centered  on  the  western 
slopes  of  the  Sierra  Nevada  Mountains  and  along  the  eastern  skirt  of  the  Sacra¬ 
mento  Valley,  between  Sacramento  and  Redding,  California. 
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the  effect  of  varying  the  NTC  site  DTED  grid  spacing  from  250 


Figure  13.  Illustrating  the  eitect  ot  varying  the  Project  WIND  site  DTED  grid  spacing  from  250  m 


For  this  research  work,  meso  7  -  micro  a  data  was  chosen  from  the  Project 
WIND  data  base.  These  data  was  collected  in  the  vicinity  of  Red  Bluff,  California 
(see  Fig.  6).  The  data  included  wind  speed  and  direction  observations  taken  from 
sensors  10  m  above  the  ground  at  6  different  sites  (S5-S10);  an  upper  air  sounding 
was  also  taken  every  two  hours  at  site  S8.  Samples  of  Project  WIND  surface  and 
upper  air  data  profiles  used  in  this  study  are  illustrated  in  Appendix  A. 

The  surface  and  upper  air  data  for  the  WSMR  and  NTC  wore  similar  in 
nature  to  the  Project  WIND  data,  except  that  they  were  not  reported  at  the 
same  frequency  as  in  Project  WIND  and  in  some  cases  had  to  be  estimated. 
The  WSMR  data  were  more  in  the  micro  a  regime,  and  the  NTC  data  were 
representative  of  meso  7  domains.  WSMR  and  NTC  meteorological  data  are 
also  described  in  more  detail  in  Appendix  A. 

The  sparse  data  scenario,  for  the  terrain  near  Pilsen,  Czechoslovakia,  did  not 
include  field  observation  data.  Instead,  theoretical  meteorological  initialization 
data  was  developed  to  simulate  remote,  unsupported  field  exorcise  situations. 
For  example,  a  single,  user  specified,  10  m  surface  wind  speed  and  direction 
was  coupled  to  an  estimated  upper  air  wind  profile  that  was  derived  from  a 
logarithmic  velocity  profile  (Arya,  1988) 

U  =  u*\  In-  (2) 

rC  *0 

where  U  is  the  mean  horizontal  wind  speed  as  a  function  of  friction  velocity, 
n*,  von  Karman’s  constant,  k,  the  roughness  parameter,  and  height  above 
a  surface,  2.  The  friction  velocity  term  is  defined  as  u'  =  \[^'  w^ere  T"  's 
surface  shear  stress,  and  p  is  the  fluid  density.  More  details  on  this  sparse  data 
init  ialization  scheme  for  Pilsen,  Czechoslovakia  are  located  in  Appendix  A. 
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2.3  The  WOCSS  Model 


The  WOCSS  model’s  approximations  of  wind  velocity  fields  combines  selected 
dynamic  and  kinematic  fluid  flow  considerations.  Kinematics  centers  on  the 
development  of  solutions  to  flow  problems  without  considering  the  actual  forces 
at  work  (Hess,  1979),  while  dynamic  meteorology  considers  a  broader  spectrum  of 
problems,  including  the  fluxes  and  conversions  of  mechanical,  as  well  as  thermal 
energy,  and  equilibrium  states  of  the  atmosphere  (Riegel,  1974). 

A  fundamental  kinematic  design  features  in  the  WOCSS  model  address  ini¬ 
tial  flow  modeling  problems  such  as  mass  consistency  and  nondivergence.  For 
example,  the  WOCSS  model  assumes  that  flow  takes  place  within  the  bounds 
of  specified  surfaces.  These  surfaces  are  somewhat  different  from  sigma  surfaces, 
which  are  solved  for  using  the  following  equation  (Pielke,  1984): 


(*  ~  *a) 


(3) 


where  <r  is  the  normalized  terrain  following  layer  being  solved  for,  s  specifies  the 
upper  flow  boundary  height,  z  is  the  layer  elevation  above  the  terrain,  and  zq 
is  the  terrain  height  above  MSL.  In  the  WOCSS  model,  wherever  a  flow  surface 
intersects  the  terrain,  wind  velocity  is  set  to  zero  (because  the  flow  field  cannot 
pass  into  the  terrain  obstacle),  and  mass  conservation  theoretically  forces  flow 
around  the  intersection  with  terrain. 

The  WCCSS  coordinate  system,  therefore,  is  a  flow-following,  not  a  terrain- 
following  system.  A  major  assumption  in  the  WOCSS  flow-following  system  is 
that  there  are  no  perpendicular  components  of  flow  to  the  bounding  flow  layer 
surfaces  (Ludwig  and  Endlich,  1988).  These  flow  bounding,  or  critical  streamline, 
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surfaces  can  be  solved  for  using  a  variation  of  the  following  equation: 


(^mM  Zq) jlow  layer 


=  V  ii-—\ 


-1/2 


(4) 


Todz0 

where  zmax  is  the  highest  altitude  limit  of  flow  in  layer  n,  za  is  the  lowest  altitude 
of  the  same  flow  layer,  V0  is  the  mean  flow  velocity  along  the  z0  surface,  g  is  the 
acceleration  due  to  gravity,  T0  is  the  mean  temperature  along  z0 ,  and  is  the 
vertical  potential  temperature  gradient.  Nondivergence  is  forced  on  all  n  flow 
layers  by  the  following  continuity  equation  approximation 


du .  dv.  _ 
dx  dy 


(5) 


where  u,  =  uAz,  and  vm  =  t>Ax.  The  u.  and  vm  terms  are  flux-related  variables, 
where  A z  is  the  average  separation  between  zQ  and  zmot,  the  u  term  is  the  wind 
velocity  component  in  the  x-axis  direction,  and  v  is  the  wind  velocity  component 
along  the  y-axis  direction. 

Initial  testing  of  the  WOCSS  model  involved  exercising  default  parameters  in 
model  subroutines.  These  sensitivity  studies  centered  on  identifying  parameter 
settings  that  optimize  the  WOCSS’s  dynamic  design  feature  (e.g.,  the  ability  to 
colve  for  lateral  forcing  effects).  Therefore,  this  work  concentrated  largely  on  the 
subroutines  that  establish  the  relationship  governing  critical  streamline  height, 


(6) 


where  kinetic  flow  energy  per  unit  mass  is  solved  for  on  the  left  side,  and  potential 
flow  energy  pei  unit  mass  is  solved  for  on  the  right  side  of  the  equation.  Uc  is  the 
flow  velocity  at  the  critical  streamline  height,  g  is  the  acceleration  due  to  gravity, 
T  is  the  mean  temperature  (in  °K)  of  the  layer  considered,  Z  is  the  height  of  the 
topographic  obstacle,  He  is  the  height  of  the  critical  streamline,  z  is  the  height 
that  air  is  lifted  to  in  the  layer,  and  |*  is  the  vertical  potential  temperature 
gradient  (Ludwig  and  Endlich,  1988). 
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2.4  The  OBWIND  Model 


OBWIND  is  a  basic  three-dimensional,  kinematic  wind  model  (Henmi  and  Tabor, 
1988),  and  does  not  consider  dynamic,  lateral  forcing  effects  in  its  solutions 
to  flow  field  approximations.  It  incorporates  a  terrain  following  <r  type  flow 
coordinate  system.  This  model  interpolates  surface  wind  initialization  data  using 
the  following  equation: 


=  (7) 

msl  m=l 

where  U (*,  j,  1)  is  the  calculated  value  at  the  grid  point  (i,j,  1),  Um  is  the  observed 
wind  at  site  m,  V>m(r)  is  and  r  is  the  grid  point  to  observation  site  distance. 
Upper  level  winds  are  interpolated  using  a  similar  equation: 


Em  En  0m(r)  ‘  0n(*)IAn,i» 
Em  EnV>m(r)-tM*) 


(8) 


where  rj>n(z)  is  the  vertical  weighting  function,  and  z  is  the  vertical  distance 
between  the  radiosonde  sounding  point  and  layer  k.  The  main  attribute  of  this 
model  is  its  mass  consistency  adjustment  scheme,  which  closely  resembles  the 
approach  taken  by  Endlich  (1967).  Endlich’s  method  is  based  on  several  linear 
partial  differential  equations.  The  equation  of  central  importance  to  Henmi’s 
model  is: 


du  dv  ...  .. 

Tx  +  Ty  = 


(9) 


where  u  is  the  west-east  wind  component,  v  is  the  south-north  wind  component, 
and  ^(i,  j)  is  the  divergence  at  each  grid  point.  Divergence  is  iteratively  reduced 
to  a  predefined  threshold,  which  imposes  mass  consistency  in  the  flow  field. 
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2.5  Case  Studies 


Since  the  primary  focus  of  this  research  work  was  centered  on  the  WOCSS  model, 
validation  of  the  OBWIND  model  was  not  addressed.  Henmi’s  OBWIND  model 
was  used  strictly  to  develop  a  relative  measure  of  the  improvement  that  the 
WOCSS  model  contributes  in  simulating  lateral  transfer  processes  in  wind  flow 
fields. 

To  establish  this  comparison,  the  WOCSS  model  and  the  OBWIND  model 
were  run  with  exactly  the  same  meteorological  and  terrain  data  for  WSMR, 
NTC,  Project  WIND,  and  Czechoslovakia),  at  three  different  Ax,  A y  grid  spac- 
ings  (100,  250,  and  1000  m).  This  comparison  of  the  WOCSS  and  OBWIND 
models  included  selected  measurements  of  execution  times  under  equivalent  ini¬ 
tialization  conditions.  The  WOCSS  model  was  also  tested  against  its  own  ability 
to  resolve  lateral  transfer  via  a  series  of  runs  with  each  meteorological  data  set 
where  observation  initialization  data  was  selectively  restricted. 

The  terrain  grid  spacing  selection  scheme  had  a  significant  influence  on  the 
selection  of  vertical  layers  examined  in  this  research  work.  The  Ax,  Ay  spacing 
scheme  was  used  to  establish  similiarity  in  the  vertical  scales  examined  in  this 
study;  three  elevated  velocity  field  surfaces  were  produced  in  each  wind  model 
run  1)  the  lowest  layer  was  coincident  with  the  10  m  surface  meteorological  wind 
velocity  observations;  2)  the  second  layer,  established  at  25  m  AGL,  was  a  factor 
2.5  x  the  first  layer’s  elevation  AGL  (e.g.,  proportional  to  the  100-m  to  the  250 
m  grid  space  resolution  change);  and  3)  the  third  layer  was  selected  at  a  factor 
of  10  x  the  first  layer,  which  meant  the  final  layer  was  100  m  AGL  (proportional 
to  the  100-1000  m  change  in  grid  space  resolution). 
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{ OBWIND1  I/O  CONFIGURATION 
WOCSS  &  GRAPHICS  SUPPORT 


Two  meteorological  conditions,  stable  and  unstable  cases,  were  selected  for 
each  terrain  site.  Criteria  used  to  select  this  data  are  discussed  in  Appendix  A. 
This  meteorological  data,  the  terrain  data  bases,  the  WOCSS  model,  and  the 
OBWIND  model  were  then  installed  in  a  structured  file  directory  system  (see 
Fig.  15)  on  a  Sun  Microsystems  SPARCstation  IPC  computer  workstation,  and 
Fig.  16  summarizes  the  plan  for  the  case  studies  that  were  carried  out  using  the 
OBWIND  and  WOCSS  models. 


Figure  15.  Mapping  of  the  file  directory  scheme  used  to  organize  the  research 
work  in  a  manner  that  prevented  the  possibility  of  overlapping  or  contaminating 
individual  sets  of  model  results. 
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2.6  Analysis  Techniques 


Tucker  and  Henmi  (1990)  reported  that  basic  mass  conserving,  diagnostic  models 
did  not  produce  significantly  better  agreement  with  actual  observation  data  than 
a  simple  interpolation  of  wind  field  data.  As  a  result,  the  OBWIND  model,  which 
is  a  basic  mass  conserving  model,  was  selected  for  use  in  this  study  as  the  base 
case  test  platform  with  which  to  compare  WOCSS  model  results;  the  validity 
of  the  WOCSS  model’s  ability  to  infer  the  effect  of  lateral  mechanical  forcing 
would  be  measured  in  terms  of  its  comparison  to  the  basic  mass  conserving 
OBWIND  model,  which  does  not  consider  the  effects  of  topography  on  flow  field 
distribution.  Therefore,  detailed,  statistical  comparisons  of  the  WOCSS  and 
OBWIND  outputs  had  to  carried  out. 

For  example,  the  standard  deviations  of  WOCSS  u  and  v  vector  components, 
<rm,  were  calculated  using 


=  {Dw„(0  -  Jr,.]’/"} 


(10) 


where  Wm(k)  is  the  value  produced  by  the  WOCSS  model  at  point  k  (Henmi 
and  Tabor,  1990).  Then  <rm  quantities  were  compared  to  the  standard  deviation 
of  OBWIND  values,  <r01  which  are  calculated  using 


where  W0(k)  is  an  OBWIND  model  value  at  point  k. 


Willmott  et  al.(1985)  and  Steyn  and  McKendry  (1988)  presented  several  sta¬ 
tistical  model  evaluation  measures  -  such  as  the  root- mean- squared  difference 
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(RMSD)  and  the  index  of  agreement  -  which  were  also  incorporated  into  the 
analysis  of  WOCSS  and  OBWIND  model  results.  The  root- mean-squared  differ¬ 
ence  can  be  expressed  as 

RMSD  =  £(Wm(i)  -  ».(())>  (12) 

where  n  is  the  number  of  gridpoints,  Wm  and  WQ  are,  respectively,  the  WOCSS 
case  and  OBWIND  case  values. 

The  index  of  agreement,  d,  is  defined  as: 

rf=l-fn(flM50),/Dl/'/|  +  iO'|)1]  (13) 

1 

where  P(  =  Wm(i)~Wl  and  OJ  =  W0(i)— Wl,  with  Wl  being  the  mean  OBWIND 
case  value. 

According  to  Steyn  and  McKendry  (1988),  this  index  of  agreement,  d,  has  a 
possible  range  between  1.0  (for  100  %  agreement)  to  0.0  (for  0  %  agreement), 
thereby  providing  an  objective,  critical  comparison  measure.  A  computer  pro¬ 
gram  that  implemented  this  agreement  index  test,  as  well  as  the  other  statistical 
tests  discussed  above,  was  developed  and  employed  to  automate  the  processing 
of  all  OBWIND-WOCSS  model  comparisons  (Appendix  B  provides  a  descriptive 
outline  of  this  computer  program). 

To  obtain  further  measures  of  statistical  significance  of  the  study  results, 
correlation  coefficients,  r,  of  WOCSS  and  OBWIND  layer  data  were  calculated 
using 


where  Wm(n)  and  Wm(n)  are  the  respective  WOCSS  and  OBWIND  data  u  and  v 
components  being  tested  over  the  number  of  data  elements  in  each  set  (containing 
n  elements),  and  W~  and  W0  are  the  respective  means  of  the  data  sets  being 
compared.  This  method  was  also  used  to  study  the  WOCSS  model’s  correlation 
with  terrain  elevation. 
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3  RESULTS 


The  WOCSS  model  was  configured  for  use  in  this  study,  and  all  OUWJNI) 
and  WOCSS  model  runs  were  executed  and  selectively  analyzed  (see  Section  2.5). 
Appendix  C  provides  a  summary  of  technical  details  related  to  OBWIND  and 
WOCSS  software  execution.  The  WSMR  site  tenain  is  fairly  symmetric  when 
compared  to  the  other  3  selected  sites,  which  facilitates  preliminary  interpretation 
of  study  results.  Therefore,  the  most  detailed  discussion  of  study  results  focuses 
on  the  WSMR  model  runs;  in  particular,  the  stable  case  comparisons  were  the 
most  interesting  in  terms  of  dynamic  effects  associated  with  terrain-atmosphere 
interaction  predicted  by  the  WOCSS  model. 

3.1  WSMR  Results 

3.1.1  WOCSS  Model  Configuration 

As  noted  earlier,  the  WOCSS  model  first  had  to  be  properly  configured  before  it 
could  be  executed,  and  an  important  detail  of  this  configuration  stage  involved 
establishing  the  shape  of  the  top  flow  layer  boundary.  The  shape  of  this  bound¬ 
ary  can  be  modified  by  a  Slope  Factor  (SLFAC)  term  in  the  WOCSS  source 
code.  Figure  17  illustrates  a  variety  of  SLFAC  configurations,  each  of  which  was 
implemented  in  the  WOCSS  model  (results  from  the  WSMR  stable  25  m  AGL 
flow  surface  are  summarized  in  Table  6). 

The  compressed  boundary  layer  slope  factor  in  Fig.  17  introduces  a  variable 
degree  of  resistance  to  vertical  displacement  as  a  function  of  terrain  elevation, 
while  the  restricted  boundary  layer  slope  factor  in  Fig.  17  severely  restricts 
vertical  flow  layer  displacement  (indicative  of  extremely  stable  conditions).  Con¬ 
versely,  the  relaxed  boundary  layer  slope  factor  in  Fig.  17  illustrates  a  condition 
in  which  there  is  no  resistance  to  flow  surface  vertical  displacement  as  a  func¬ 
tion  of  terrain  elevation  (more  characteristic  of  neutral- unstable  boundary  layer 
conditions). 
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Note  in  Fig,  17  that  flow  surfaces  below  the  boundary  top  can  intersect 
terrain  (e.g.,  layer  n),  where  kinetic  flow  energy  is  predicted  to  be  insufficient 
in  overcoming  local  terrain  impediments  along  the  z  axis.  Figure  18  illustrates 
the  flow  restricted  space  that  results  when  a  flow  layer  surface  intersects  the 
underlying  terrain.  The  WOCSS  model  assigns  zero  wind  speeds  to  all  grid  points 
in  flow  restricted  spaces,  which  then  requires  the  WOCSS  mass  conservation 
algorithm  to  force  flow  mass,  on  each  modeled  flow  layer,  around  intersections 
with  flow  restricted  spaces.  WOCSS  wind  speeds  are  set  to  zero  in  flow  restricted 
spaces  because  flow  in  the  vicinity  of  flow  restricted  spaces  is  predicted  to  have 
insufficient  kinetic  energy  to  overcome  the  terrain  obstacle  in  question. 

In  Table  6,  the  difference  in  flow  restricted  space  between  the  relaxed  (layer 
n+1  in  Fig.  17)  and  compressed  (layer  n+2  in  Fig.  17)  slope  factor  in  WOCSS 
model  25  m  layer  runs  is  equivalent  to  a  2.05  x  2.05  km2  area,  and  the  difference 
between  the  compressed  and  restricted  slope  (layer  n+3  in  Fig.  17)  model  runs  is 
equivalent  to  a  2.96  x  2.96  km2  area.  The  reduction  in  predicted  flow  restricted 
space  at  25  m  between  compressed  slope  runs,  after  being  re-initialized  with  wind 
speed  data  increased  by  a  factor  of  2x,  is  equivalent  to  a  3.07  x  3.07  km 2  area. 

The  differences  in  flow  restricted  space  described  in  the  previous  paragraph 
are  all  noteworthy,  especially  in  terms  of  urban  scale  problems,  and  the  differences 
are  reasonably  consistent  with  the  conditions  that  produced  them.  For  example, 
as  the  slope  of  the  boundary  layer  top  became  more  restricted,  the  amount  of 
restricted  flow  space  decreased,  because  Bernoulli  speed-up  effects  in  the  more 
compressed  flow  layers  produced  increases  in  available  kinetic  energy,  improving 
local  flow  potential  to  overcome  terrain  obstacles;  this  also  explains  the  resulting 
decrease  in  flow  restricted  space  when  the  magnitude  of  wind  speed  used  to 
initialize  the  WOCSS  model  was  doubled. 
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n  -  25  m  AGL  in  Fig.  17.  The  perimeter  of  the  flow  restricted  space  actually 
intersects  the  underlying  terrain  surface.  The  WOCSS  model  does  not  allow 
mass  transport  through  these  flow  restricted  spaces. 
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Table  6.  Cumulative  area  of  flow  restricted  spaces  (in  m3  and  %  of 
total  plane  area  (410,062,500  m3)),  produced  by  the  WOCSS  model, 
at  25  m  above  ground  level  under  stable  conditions  using  the  WSMR 
site  terrain  (see  Figs.  19-21).  Variations  in  the  slope  along  the  top  of 
the  modeled  boundary  layer  produced  significant  changes  in  the  area 
of  flow  restricted  spaces;  changing  the  speed  of  initialization  data  also 
produced  significant  variations  in  flow  restricted  spaces. 


Slope  Factor 

Stable 

Wind  speeds  doubled 

Relaxed 

31,437,500  m3  (7.7%) 

Compressed 

27,250,000  m3  (6.7%) 

17,812,500  m3  (4.3%) 

22,687,500  m3  (5.5%) 

Figures  19-21  illustrate  the  relative  areas  and  changes  in  the  25-m  flow  re¬ 
stricted  spaces  under  stable  conditions  noted  in  Table  6.  Under  the  conditions 
examined  in  Table  6,  unstable  WOCSS  model  runs  did  not  produce  any  flow 
restricted  spaces  at  10,  25,  or  100  m  above  ground  level,  because  the  increased 
average  speeds  in  the  unstable  case,  only  2.45  m/s  greater  than  in  the  stable  case 
at  10  m  AGL  (see  Appendix  A),  provided  sufficient  kinetic  energy  to  selected 
flow  surfaces  to  overcome  all  intervening  WSMR  site  terrain  obstables.  This  re¬ 
sult,  in  conjunction  with  the  stable  case  results,  illustrates  the  sensitivity  of  the 
WOCSS  model  to  atmospheric  initialization  conditions.  Because  this  sensitivity 
is  indicative  of  processes  actually  observed  in  nature,  this  WOCSS  model  respon¬ 
siveness  is  a  significant  feature  that  will  contribute  positively  to  WOCSS  model 
application  potential. 
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Figure  20.  Same  as  Fig.  19,  except  that  a  compressed  top  boundary  slope  factor  (layer  n+2  in  Fig.  17) 
was  used  to  predict  the  illustrated  flow-restricted  regions. 
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Figure  21 .  Same  as  Fig.  20.  except  that  all  surface  and  vertical  wind  initialization  speeds  were 
increased  by  a  factor  of  2X  the  original  observed  speed  magnitudes  listed  in  Appendix  A. 


For  the  purposes  of  this  study  a  compressed  slope  factor  (SLFAC=0.5)  was 
selected  for  use  in  the  stable  case  model  runs,  and  a  relaxed,  boundary  layer 
slope  factor  (SLFAC=1.0)  was  used  in  unstable  model  runs.  Under  stable  con¬ 
ditions  the  atmosphere  is  more  resistant  to  vertical  displacement,  and  this  effect 
is  modeled  along  flow  layers  over  elevated  surfaces  when  the  boundary  layer  top 
is  compressed. 

The  compression  referred  to  in  the  previous  paragraph  results  in  selective 
constriction  of  flow  layer  thickness  in  the  vicinity  of  elevated  surfaces,  and  thereby 
simulates  Bernoulli  speed-up  effects  over  relative  elevation  maximums  (layer  n+2 
in  Fig.  17).  Conversely,  a  relaxed  boundary  layer  slope  factor  (layer  n+1  in  Fig. 
17)  was  used  in  unstable  WOCSS  model  runs  to  simulate  the  effect  of  reduced 
resistance  to  vertical  displacement,  which  is  consistent  with  unstable  conditions 
observed  in  nature. 

3.1.2  OB  WIND- WOCSS  Statistical  Comparisons 

The  minimum,  mean,  maximum,  and  ±  1.0  standard  deviation  of  each  stable 
case  WSMR  model  run  axe  listed  in  Table  7  (u  components)  and  Table  8  (v 
components).  Summary  statistics  (e.g.,  the  root  mean  square  difference,  agree¬ 
ment  index,  and  correlation  coefficient)  are  then  listed  in  Table  9  (u  components) 
and  Table  10  (v  components). 

The  u  component  represents  approximately  60%  of  the  mean  flow  vector  over 
the  WSMR  site  stable  case  model  runs.  Therefore,  significant  OBWIND-WOCSS 
comparison  differences  are  more  apparent  in  the  u  component  statistics.  Note 
the  0.0  WOCSS  minimum  u  components  in  Table  7  at  10  and  25  m.  These  0.0 
minimums  are  associated  with  the  flow  restricted  spaces  discussed  earlier  (also 
see  Figs.  19-21).  Table  8  also  lists  0.0  WOCSS  v  component  minimums  at  10 
and  25  m. 
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Grid  spacing  selection  does  not  appear  to  significant!}'  afFect  minimum  u 
components.  For  example,  75%  of  all  OBWIND  and  WOCSS  non-zero  min¬ 
imum  u  components  consistently  decreased  an  average  of  only  3.12%  as  grid 
resolution  increased  from  1000-100  m.  However,  large  differences  occurred  be¬ 
tween  WOCSS  and  OBWIND  non-zero  tninimums  at  100m  (288%  change  be¬ 
tween  100  m  grid  OBWIND-WOCSS  minimum  u  components,  dropping  only  to 
a  260%  change  at  the  1000  m  grid  stable  runs,  with  OBWIND  producing  smaller 
non-zero  u  and  v  minimums). 

Mean  10  m  OBWIND  u  component  magnitudes  were  196%,  192%,  and  157% 
greater  than  WOCSS  means  at  100,  250,  and  1000  m  grid  spacing,  respectively, 
but  this  completely  reversed  in  the  25  and  100  m  level  outputs  (for  example,  at 
25  m  WOCSS  means  were  114%,  115%,  and  124%  greater  than  OBWIND  means 
over  the  same,  respective  grid  spacing  selection).  At  the  10  m  level  WOCSS 
means  increased  by  15.6%  and  18.4%,  respectively,  as  grid  spacing  resolution 
increased  from  1000-250,  and  1000-100  m. 

Maximum  10  m  OBWIND  u  components  were  also  greater  than  WOCSS 
maximums  by  50%,  42%,  and  11%,  respectively,  at  100,  250,  and  1000  m  grid 
spacing,  and  as  in  the  case  of  mean  component  s,  WOCSS  maximums  were  greater 
than  OBWIND  maximums  at  the  25  and  100  m  levels  (by  179%,  205%,  and  250% 
at  100,  250,  and  1000  m  grid  spacing,  respectively). 

Maximum  OBWIND  u  and  v  components  were  also  always  higher  when  the 
100  m  grid  was  employed;  conversely,  maximum  WOCSS  u  components  were 
always  higher  when  the  1000  m  grid  was  employed  (v  component  maximums  be¬ 
haved  erratically.  These  differences  are  attributed  to  variations  in  OBWIND  and 
WOCSS  interpolation  algorithms  and  how  these  algorithms  assign  wind  velocity 
values  to  individual  grid  points  at  different  grid  spacing  scales. 
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Table  7.  Basic  u  component  windfield  statistics  for  stable  WSMR 
site  OBWIND  (u0)  and  WOCSS  (uw)  windfield  outputs  (min=Minimum, 
avg=Mean,  max=Maximum,  and  a  =±  1  Standard  Deviation). 


Model  Source  Statistic  100m  grid  f  250m  grid 

1000m  grid 

u0 

min  10m 

■MEW 

■»y— 

-0.716014 

uw 

min  10m 

0.000000 

0.000000 

0.000000 

U0 

min  25m 

-0.456276 

uw 

min  25m 

0.000000 

0.000000 

0.000000 

U„ 

min  100m 

-0.462111 

■TiTO-W 

-0.487063 

uw 

min  100m 

-1.796830 

-1.743190 

-1.752220 

u  o 

avg  10m 

•0.984218 

-0.969837 

-0.902576 

Uw 

avg  10m 

-0.331429 

^■lcLiinx^i 

Uo 

avg  25m 

-0.562656 

-0.554448 

Uw 

avg  25m 

-1.201638 

-1.195562 

-1.161671 

u„ 

avg  100m 

wmibwm 

mxmmm 

-0.528225 

Uw 

avg  100m 

-1.898141 

-1.878264 

-1.793024 

Uo 

BufTSTiHiV 

-1.675720 

-1.602960 

-1.400490 

Uw 

max  10m 

-1.114420 

-1.126640 

-1.258820 

Uo 

max  25m 

TiKffliTiEy 

-0.746207 

-0.708540 

Uw 

max  25m 

-1.864310 

-1.851610 

-1.922880 

Uo 

max  100m 

mswnmm 

-0.711103 

Uw 

max  100m 

-2.202610 

-2.278180 

-2.486910 

Uo 

a  10m 

0.084747 

0.087054 

0.119891 

Uw 

a  10m 

0.458096 

0.462548 

0.486732 

Uo 

a  25m 

0.041628 

0.038673 

0.063885 

Uw 

a  25m 

0.332401 

0.337656 

0.391335 

Uo 

<r  100m 

0.038507 

0.036918 

0.063807 

Uw 

o  100m 

0.024130 

0.053976 

0.184966 
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Table  8.  Basic  v  component  windfield  statistics  for  stable  WSMR 
site  OBWIND  (v„)  and  WOCSS  (vtt)  windfield  outputs  (min=Minimum, 
avg=Mean,  max=Maximum,  and  a  =±  1  Standard  Deviation). 


Model  Source  Statistic  100m  grid  250m  grid 

1000m  grid 

V* 

min  10m 

■TiEHfrWl 

vw 

■mWOTThW 

0.000000 

0.000000 

0.000000 

Vo 

mmmm 

vw 

min  25m 

0.000000 

0.000000 

0.000000 

V0 

min  100m 

-0.250070 

-0.287582 

-0.330255 

Vw 

min  100m 

-1.264500 

-1.231830 

-1.275420 

v„ 

-0.625345 

-0.615236 

-0.567920 

Vw 

avg  10m 

-0.238849 

mivxrmw 

-0.243642 

v0 

avg  25m 

-0.379551 

-0.346868 

Vw 

avg  25m 

mmuzm 

-0.855768 

V0 

avg  100m 

-0.373936 

m*mmM 

-0.341767 

Vw 

avg  100m 

-1.377465 

-1.359816 

-1.287632 

Vo 

max  10m 

-0.935505 

-0.932253 

V» 

max  10m 

-0.845871 

Vo 

max  25m 

-0.514950 

-0.485605 

-0.445452 

Vw 

max  25m 

-1.626910 

-1.493900 

-1.543850 

V0 

max  100m 

-0.502560 

-0.476896 

-0.438018 

Vw 

max  100m 

-1.687400 

-1.723340 

-1.679680 

Vo 

a  10m 

0.115050 

0.114872 

0.122096 

Vw 

a  10m 

0.330129 

0.329335 

0.336934 

Vo 

a  25m 

0.053658 

0.045551 

0.042943 

Vw 

o  25m 

0.240193 

0.243583 

0.283064 

Vo 

a  100m 

0.050773 

0.043380 

0.041844 

Vw 

a  100m 

0.018221 

0.039398 

0.130902 
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Table  9.  Summary  u  component  windfield  statistics  for  WSMR  site 
10, 25,  and  100  m  stable  case  OBWIND  (u0)  and  WOCSS  (uw)  wind- 
field  outputs  (rms=Root  Mean  Square  Difference,  agi=Agreement 
Index,  and  r=Correlation  Coefficient). 


Component  Sources 

Statistic  100m  grid  |  250m  grid  1000m  grid 

u0*uw 

rms  10m 

0.809818 

0.806546 

0.781163 

u0-uw 

rms  25m 

0.730285 

0.741253 

0.812019 

u0-uw 

rms  100m 

1.337306 

1.345582 

1.395853 

u0-uw 

agi  10m 

0.092185 

0.074378 

0.056239 

U0-Uw 

agi  25m 

0.085009 

0.104984 

0.209522 

U0-Uw 

agi  100m 

0.048872 

0.060108 

0.135892 

Uo-Uu, 

r  10m 

-0.068655 

-0.059891 

0.051799 

U„-U  w 

r  25m 

-0.113005 

0.068513 

0.346386 

u0-u«, 

r  100m 

0.009543 

0.223450 

0.742598 

Table  10.  Summary  v  component  windfield  statistics  for  WSMR  site 
10,  25,  and  100  m  stable  case  OBWIND  (vn)  and  WOCSS  (vw)  wind- 
field  outputs  (rms=Root  Mean  Square  Difference,  agi=Agreement 
Index,  and  r=Correlation  Coefficient). 


Component  Sources  1  Statistic 


rms  10m 


25m 


rms  100m 


250m  grid  1000m  grid 


r  10m 


r  25m 


r  100m 


0.506293 


0.551708 


1.014435 


0.268167 


0.147221 


0.083071 


0.242422 


0.504318 


0.555286 


1.017097 


0.248396 


.141707 


0.242713 


0.069937 


0.204643 


0.487292 


0.588299 


1.041666 


0.201858 


0.195270 


0.126815 


0.276910 


0.174258 


0.785983 
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Table  11.  Summary  u  component  windfield  statistics  for  WSMR  site 
10,  25,  and  100  m  unstable  case  OBWIND  (u0)  and  WOCSS  (uw) 
windfield  outputs  (rms=Root  Mean  Square  Difference,  agi= Agreement 
Index,  and  r=Correlation  Coefficient). 


1  Component  Sources  Statistic  100m  grid  250m  grid  1000m  grid 

u0-uw 

|  rms  10m 

2.094479 

2.138106 

u0-u„ 

[  rms  25m 

0.674350 

0.613055 

0.537859 

Uo-U,,, 

|  rms  100m 

0.245543 

0.190512 

0.183246 

u0-uw 

agi  10m 

0.207519 

0.135198 

0.111055 

Uo-Uu, 

agi  25m 

0.273099 

0.195771 

0.126571 

u0-uw 

agi  100m 

0.412048 

0.312226 

0.487617 

Uo-Uu, 

r  10m 

0.087668 

Mimmm 

0.307092 

UoUu, 

r  25m 

-0.026439 

0.135907 

0.471575 

Uo'Uu, 

r  100m 

0.012076 

0.224244 

0.613856 

Table  12.  Summary  v  component  windfield  statistics  for  WSMR  site 
10,  25,  and  100  m  unstable  case  OBWIND  (v0)  and  WOCSS  (vj 
windfield  outputs  (rms=Root  Mean  Square  Difference,  agi= Agreement 
Index,  and  r=Correlation  Coefficient). 


Component  Sources  Statistic  100m  grid  )  250m  grid  |  1000m  grid 

Vo-Vw 

rms  10m 

0.759228 

1.230143 

0.785934 

v0-vw 

rms  25m 

1.266384 

1.326321 

1.216271 

V„-Vw 

rms  100m 

2.597788 

2.740088 

2.606363 

Vo-Vw 

agi  10m 

0.134758 

0.328341 

0.212271 

Vo-Vw 

agi  25m 

0.067452 

0.258935 

0.351719 

Vo-Vw 

KfflTiTij'y 

0.033449 

0.137587 

0.196585 

Vo-Vu, 

r  10m 

0.058383 

0.128353 

0.403180 

V0-Vw 

r  25m 

0.406679 

0.173538 

0.943587 

Vo-Vw 

r  100m 

0.435636 

0.183103 

0.955592 
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In  general,  the  WOCSS  model  v  components  (Table  8)  behaved  in  a  similar 
manner  to  the  u  components,  producing  significantly  greater  winds  than  014- 
WIND  at  25  and  100  m  levels  (outside  of  flow  restrict'd  spaces);  and  the  WOCSS 
ouput  winds  had  significantly  greater  variance  in  distributions  at  10  and  25  m 
(see  the  standard  deviation  data  in  Tables  7  and  8).  Variance  in  mean  WOCSS 
standard  deviations  also  increased  as  the  grid  spacing  increased  (by  17.7%  from 
100-1000  m  and  15.9%  from  250-1000  m  grid  spacing).  The  increased  variance 
at  10  and  25  m  is  exaggerated  because  of  the  effect  of  flow  restricted  space  in  the 
WOCSS  model  outputs  at  10  and  25  m  in  the  stable  case  scenario. 

Comparison  of  the  summary  statistics  in  Tables  9-10  (stable  case)  and  Tables 
11-12  (unstable  case)  indicate  very  low  agreement  indices  and  correlation  between 
OBWIND  and  WOCSS  model  results  in  89%  of  all  cases  (e.g.,  average  u  and  v 
component  agreement  indices  were,  respectively,  0.096353  and  0.198378  in  the 
stable  case,  and  0.25123-1  and  0.191233  in  the  unstable  case);  unstable  agreement 
was  50%  better  than  stable  case  agreement. 

Only  the  100  m  level  u  and  v  components  with  1000  m  grid  spacing  had  consis¬ 
tently  good  correlation  (r  >  0.5),  but  this  high  correlation  case  was  not  associated 
with  the  largest  agreement  index.  Therefore,  agreement  indices  and  correlation 
coefficients  do  not  appear  to  produce  proportionally  equivalent  results.  However, 
agreement  and  correlation  generally  deteriorated  as  grid  resolution  increased. 

Tables  9-12  imply  that  agreement  indexes  and  correlation  coefficients  arc  gen¬ 
erally  higher  for  the  unstable  case,  where  more  kinetic  energy  was  available  in  the 
initialization  data  (winds  at  10  m  in  the  unstable  case  were  approximately  10-1% 
greater  in  magnitude  than  in  the  stable  case  (see  Appendix  A  for  a  listing  of  the 
stable  and  unstable  initialization  data)).  Because  of  the  availability  of  sufficient 
kinetic  energy  in  the  unstable  case,  no  ilow  restricted  space  was  predicted  by  the 
WOCSS  model  at  10,  25,  and  100  m  AGL,  which  largely  explains  the  improved 
agreement  and  correlation  with  OBWIND  model  results  compared  to  the  WSMK 
stable  case  comparisons. 
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3.1.3  OBWIND-WOCSS  Graphical  Comparisons 

Results  presented  up  to  this  point  have  been  in  tabular  form,  and  significant 
results  have  been  referred  to  only  through  textual  description  (primarily  to  econ¬ 
omize  on  the  total  number  of  graphical  plots  that  would  be  necessary  to  display 
all  the  previously  discussed  relationships  among  the  WSMR  results).  However, 
graphical  presentation  of  selected  model  output  fields  becomes  necessary  at  this 
point  to  illustrate  the  value  of  spatial  analysis  of  model  results,  which  in  many 
cases  provides  more  useful  information  about  model  performance  than  any  of  the 
previously  discussed  statistics. 

For  example,  Fig.  22  is  a  wind  vector  field  plot  of  the  OBWIND  10  m  stable 
output  at  250  m  grid  spacing  (81  x  81  grid  points).  The  density  of  this  data 
output  nearly  saturates  the  plot  area  and  was  therefore  used  as  the  useful  upper 
threshold  in  producing  plots  of  model  outputs  (e.g.,  vector  plots  of  100  m  grid 
model  outputs  were  illegible  in  many  cases). 

Figure  23  is  a  qualitative  graphical  representation  of  stable  case  OBWIND 
wind  speeds,  (\/ua  +  u2),  associated  with  the  10  m  layer  wind  vectors  plotted  in 
Fig.  22.  Figures  24  and  25  are  the  10  m  WOCSS  plots  corresponding  to  the 
OBWIND  plots  in  Figs.  22  and  23  respectively.  Notice  the  zeroing  effect  of  the 
flow  restricted  space  predicted  by  the  WOCSS  model  in  the  center  of  Fig.  24. 
This  effect  is  more  obvious  in  Fig.  25,  where  speed  is  plotted  as  a  surface  over 
the  WSMR  terrain. 

Wind  vector  differences,  A  =  W’wocss  -  Vobwind »  were  also  calculated  (see 
Fig.  26).  Negative  difference  vectors  are  indicated  in  the  center  of  the  wind- 
field,  which  coincide  with  the  WOCSS  flow  restricted  space  over  higher  terrain 
elevations  (see  Figs.  19-21).  Outside  of  the  flow  restricted  spaces  direction  agree¬ 
ment  with  OBWIND  results  are  good,  but  WOCSS  wind  speeds  are  generally 
higher  than  OBWIND  wind  speeds  over  the  entire  unrestricted  flow  space  (most 
noticable  at  25  and  100  m  AGL).  This  is  attributed  to  WOCSS  mass  balancing 
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OBWIND  1 0M  STABLE  250M  GRID 
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Figure  22.  OBWIND  10  m  stable  case  wind  vector  field  plot  for  the  WSMR  site  using  250  X  250  square  m  grid 
cells  (81  X  81 ).  WSMR  site  terrain  elevation  contours  are  also  included  (north  is  located  along  the  top,  labeled 
edge  of  the  plot).  Magnitude  of  minimum,  average,  and  maximum  speed  vectors  are,  respectively,  0.917, 1.1 54, 
and  1.638  m/s  (vector  plot  size  is  proportional  to  the  speed  associated  with  each  vector). 


spEEO 
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Figure  24.  WOCSS  10  m  stable  case  wind  vector  field  plot  for  the  WSMR  site  using  250  X  250  square  m 
grid  cells  (81  X  81).  WSMR  site  terrain  elevation  contours  are  also  included  (north  is  located  at  the  top  labeled 
edge  of  the  plot).  Minimum,  average,  and  maximum  vector  speed  magnitudes  are,  respectively,  0.000, 0.408, 
and  1 .390  m/s  (vector  plot  size  is  proportional  to  the  speed  associated  with  each  vector). 
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WOCSS-OBWIND  1 0  m  stable  case  wind  vector  difference  field  plot  for  the  WSMR  site  usina 
250  X  250  square  m  grid  cells  (81  X  81).  WSMR  site  terrain  elevation  contours  are  also  included  (north  is 
located  at  the  top  labeled  edge  of  the  plot).  Magnitude  of  the  average  vector  difference  «  0  790  m/s  and 
magnitude  of  the  maximum  vector  difference  -  1 .63 8  m/s. 


adjustments  over  the  remaining  flow  space,  which  provides  gross  adjustment  for 
the  effect  of  flow  restricted  spaces  on  the  windfield  mass  consistency  via  speed 
adjustments  to  all  free  flow  grid  volumes. 

The  result  described  in  the  previous  paragraph  is  consistent  with  data  listed 
in  Tables  7  and  8;  it  is  also  reflected  in  the  behavior  of  WOCSS  wind  speed- 
terrain  correlation.  For  example,  using  100  m  grid  spacing,  at  10  and  25  m  AGL 
the  WOCSS  wind  speed  correlation  with  underlying  WSMR  terrain  under  stable 
conditions  is  very  poor  (r=-0. 034921  and  -0.068222,  respectively).  However,  at 
100  m  AGL  under  the  same  circumstances  correlation  is  very  strongly  inversely 
related  (r=-0.779617). 

At  10  and  25  m  AGL  flow  restricted  space  diminishes  correlation  of  the 
WOCSS  wind  fields  with  terrain,  but  at  100  m  AGL  (where  there  was  no  pre¬ 
dicted  flow  restricted  space)  the  generally  monotonic  winds  are  more  inversely 
related  to  the  lower  elevations  than  the  higher  elevations  (primarily  because  a 
majority  of  the  WSMR  terrain  was  composed  of  a  low  basin  plain,  the  inverse 
relationship  of  WOCSS  100  m  winds  with  terrain  was  strong).  This  inverse  cor¬ 
relation  would  not  have  been  as  significant  at  100  m  AGL  if  the  WOCSS  model 
predicted  localized  Bernoulli  speed  up  effects  in  the  vicinity  of  higher  terrain. 

The  wind  vector  difference  plots  illustrate  that  the  WOCSS  mass  balancing 
algorithm  does  not  calculate  dyjdx  or  dx/dy  horizontal  wind  direction  changes 
that  would  result  in  nature  near  flow  restricted  spaces,  and  the  WOCSS  model 
also  does  not  calculate  local  horizontal  Bernoulli  speed  up  effects,  unless  these 
effects  were  incorporated  in  the  initialization  wind  data.  These  are  all  significant 
limitations  to  the  WOCSS  model’s  potential  to  produce  valid,  realistic  windfield 
outputs.  However,  these  limitations  could  be  compensated  for  to  varying  degrees 
by  selectively  coupling  WOCSS  output  with  other  models  (which  will  be  discussed 
in  Section  4). 
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Figure  27  is  a  plot  of  25  m  OB  WIND  stable  case  WSMR  data  at  1000  m 
x  1000  m  grid  spacing  and  is  followed  by  a  plot  of  the  same  data  in  Fig.  28, 
where  grid  spacing  has  been  set  to  250  m  x  250  m;  Figures  27-28  do  not  show 
significant  differences  in  grid  spacing  effect  because  of  the  monotonic  behavior  of 
the  general  wind  field  resulting  from  initializing  the  model  with  only  one  vertical 
wind  profile.  However,  when  the  OBWIND  speed  surface  plots  are  compared 
(Figs.  29  and  30)  using  similar  grid  spacing  selections,  subtle  speed  maximums 
are  apparent  in  the  250  m  grid  plot  (Fig.  30)  that  are  missed  in  the  1000  m  grid 
plot  (Fig.  29). 

Figures  31-34  illustrate  the  same  sequence  of  plots  using  25  m  WOCSS  stable 
case  outputs,  which  compare  to  Figs.  27-30  respectively.  Note  that  flow  restricted 
space  at  25  m  AGL  is  significantly  smaller  than  that  at  10  m  AGL  (refer  to  Figs. 
24  and  32);  WOCSS  flow  is  less  influenced  by  terrain  effects  as  flow  surface 
elevation  above  ground  level  increases,  as  observed  in  nature. 

The  25  m  level  vector  differences  are  also  contrasted  using  250  and  1000  m  grid 
plots  in  Figs.  35  and  36,  respectively.  Comparison  of  the  center  region  differences 
in  the  vicinity  of  the  highest  terrain  illustrate  a  smoothing  effect  in  the  1000  m 
grid  plot  (Fig.  35);  the  250  m  gridded  model  output  would  be  significantly  more 
accurate  in  performing  high  resolution  flow  trajectory  operations  (for  example, 
in  urban  scale  applications). 

Figures  37-39  compare  WSMR  site  100  m  stable  case  OBWIND  and  WOCSS 
model  outputs  using  250  m  grid  spacing.  In  this  comparison  the  WOCSS  model 
predicts  no  flow  restricted  spaces  at  100  m  AGL  (Fig.  38),  and  WOCSS  winds 
are  all  greater  than  those  predicted  by  the  OBWIND  model  (Fig.  39).  These 
results  are  reflected  in  Tables  12  and  13  but  are  more  obvious  in  graphical  spatial 
presentations  (e.g.,  Figs.  37-39). 
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Figure  28.  OBWIND  25  m  stable  case  wind  vector  field  plot  for  the  WSMR  site  using  250  X  250  square  m 
grid  cells  (81  X  81 ).  WSMR  site  terrain  elevation  contours  are  also  included  (north  is  located  at  the  top  labeled 
edge  of  the  plotted  regions).  Magnitude  of  the  average  speed  vector «  0.670  m/s,  and  magnitude  of  the 
maximum  speed  vector  -  0.821  m/s. 
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Figure  31 .  WOCSS  25  m  stable  case  wind  vector  plot  for  the  WSMR  site  using  1000  X  1000  square  m 
grid  ceils  (21  X  21 ).  WSMR  site  terrain  elevation  contours  are  also  included  (north  is  located  at  the  top 
labeled  edge  of  the  plot).  Magnitude  of  the  average  speed  vector  -  1 .423  m/s,  and  magnitude  of  the 
maximum  speed  vector  «  2.139  m/s. 
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Figure  32.  WOCSS  25  m  stable  case  wind  vector  field  plot  for  the  WSMR  site  using  250  X  250  square  m 
grid  ceils  (61  X  81 ).  WSMR  site  tenain  elevation  contours  are  also  included  (north  is  located  at  the  top  labeled 
edge  of  the  plot).  Magnitude  of  average  speed  vector  - 1 .472  m/s,  and  magnitude  of  the  maximum  speed 
vector  -  2.070  m/s. 
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WOCSS— OBWIND  VECTOR  DIFFERENCE  25M  STABLE  1000M  GRID 
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Figure  35.  WOCSS-OBWIND  25  m  stable  case  wind  vector  diflerence  field  plot  for  the  WSMR  site  using 
1000  X 1000  square  m  grid  cells  (21  X  21 ).  WSMR  site  terrain  elevation  contours  are  also  included  (north 
is  located  at  the  top  labeled  edge  of  the  plot).  Magnitude  of  the  average  vector  differences  «  0.892  m/s, 
and  magnitude  of  the  maximum  vector  diffemce  « 1 .454  m/s. 
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Figure  36.  WOCSS-OBWIND  25  m  stable  case  wind  vector  differences  plot  for  the  WSMR  site  using 
250  X  250  square  m  grid  cells  (81  X  81).  WSMR  site  terrain  elevation  contours  are  also  included  (north 
is  located  at  the  top  labeled  edge  of  the  plot).  Magnitude  of  average  vector  differences  »  0.892  m/s, 
and  magnitude  of  the  maximum  vector  difference  - 1 .328  m/s. 


OBWIND  mOM  STABLE  25QM  GRID 


WOCSS  100M  STABLE  250M  GRID 


m 

>>X*N>\V 


\\VA^V§VW\ 

»■! 


liiiililli 


.2>e-5' 

u.  o>  ® 


70 


lure  38.  WOCSS  100  m  stable  case  wind  vector  field  plot  for  the  WSMR  site  using  250  X  250  square  m 
i  cells  (81  X  81 ).  WSMR  site  terrain  elevation  contours  are  also  included  (north  is  located  at  the  top  lab 
3®  of  the  plot).  Magnitude  of  average  speed  vector  =  2.319,  and  the  maximum  speed  vector  =  2.745  m/: 
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Figure  39.  WOCSS-OBWIND  100  m  stable  case  wind  vector  difference  field  plot  for  the  WSMR  site  using 
250  X  250  square  m  grid  cells  (81  X  81).  WSMR  site  terrain  elevation  contours  are  also  included  (north  is 
located  at  the  top  labeled  edge  of  the  plot).  Magnitude  of  average  vector  differences  «  1 .642  m/s,  and  the 
magnitude  of  the  maximum  vector  difference  «  2.077  m/s. 


Figure  40  illustrates  the  completely  non  zero  WOCSS  unstable  10  in  layer 
wind  vector  field.  All  unstable  WOCSS  layer  outputs  were  completely  non-zero, 
resulting  from  the  increased  availability  of  kinetic  energy  in  each  layer.  Because 
of  the  completely  free  flow  nature  of  the  unstable  case,  graphical  plots  of  results 
are  generally  very  monotonic  and,  therefore,  provide  little  additional  analysis 
value  over  the  summary  statistics  (Tables  1C  and  17). 

3.1.4  Initialization  Data  Correlation  Test 

The  final  test  involved  comparing  WOCSS  and  OBW1NI)  I  aver  outputs  with 
varying  degrees  of  surface  initialization  data  (e.g.,  with  nm,ir.  ii„„u  -  1,  i - 
2,  ...  surface  observations).  Generally  this  testing  produced  insignificant  differ¬ 
ences  in  the  WOCSS  model  outputs.  For  example,  Table  13  indicates  stable  case 
WOCSS  u  component  agreement  index  values  were  never  below  1.0;  WOCSS 
correlation  coefficients  (Table  14)  also  never  fell  below  1.0  in  all  10  m  stable 
cases  (the  level  at  which  all  surface  observations  were  actually  acquired). 

The  high  WOCSS  and  OBWIND  agreement  index  f  A  correlation  coefficient 
values  produced  in  this  testing  served  to  illustrate  the  upper  range  ielial>ilify 
of  both  the  agreement  index  and  correlation  coefficient  algorithms  listed  in  Ap¬ 
pendix  B  (e.g.,  agreement  index  and  correlation  coefficient  values  must  he  <  1 .()). 
However,  these  results  also  imply  that  the  version  of  the  WOCSS  model  being 
evaluated  does  not  consider  supplemental  surface  based  observation  data  dining 
wind  field  interpolation  operations;  apparently  only  the  vertical  wind  profile  data 
is  during  these  WOCSS  operations,  which  also  largely  explains  why  the  WOCSS 
model  run-time  is  significantly  less  than  for  the  OBWIND  mod<  •1  (see  Appendix 

C). 
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The  WOCSS  model’s  emphasis  on  weighting  upper  air  data  is  also  large  ly 
responsible  for  the  greater  root  mean  square  differences  in  the  unstable  WSMH 
cross-model  comparisons  listed  in  Tahles  9  12  (e.g.,  in  the  unstable  case,  the 
WOCSS  model  did  not  consider  the  increased  low  level  variability  in  surface 
conditions  while  the  OBWIND  model  did  (see  Appendix  A  to  cont  rast  the  WSM  H 
stable  and  unstable  surface  observation  variability)). 

Variation  in  the  OBWIND  10  m  agreement  index  and  correlation  coefficient 
values  illustrates  several  points.  The  most  important  point  is  that  the  OBWIND 
model  is  more  responsive  to  surface  initialization  data  than  the  WOCSS  model. 
The  second  point  involves  the  last  2  OBWIND  values  in  Tables  1.1  and  H;  in 
these  2  cases,  agreement,  index  and  correlat  ion  coefficients  are  not  i  onsistenl  wit  li 
each  other.  This  result  makes  it  clear  that  dependency  on  a  single  statistical 
evaluation  would  probably  be:  misleading  in  certain  cases. 

The  basic  performance  characteristics  of  both  the  OBWIND  and  WOCSS 
models  are  evident  from  review  of  WSMR  results.  Therefoir,  the-  remaining 
results  will  be  expanded  upon  only  when  necessary  and  will  primarily  locus  on 
performance  characteristics  of  the  WOCSS  model,  emphasizing  details  regarding 
the  site  portability  aspects  of  the  WOCSS  model  with  respect  to  variation  of 
terrain  complexity. 
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Table  13.  Summary  of  OBWIND  and  WOCSS  model  u  component 
agreement  index  value  (agi)  comparisons  for  the  WSMR  site  10  m 
stable  case  un  windfieid  outputs  where  n  is  the  number  of  surface 
observations  used  to  initialize  model  runs  with  a  250  m  grid.  Note 
that  the  results  presented  here  were  produced  by  testing  each  model 
against  itself. 


Component  Sources  Statistic  OBWIND  WOCSS 

u$-u5 

agi 

1.000000 

1.000000 

U4-U5 

agi 

0.437999 

U3-U5 

agi 

0.458919 

1.000000 

U2-U5 

agi 

0.509499 

1.000000 

Ui-U5 

agi 

0.079589 

1.000000 

Table  14.  Summary  of  u  component  correlation  coefficient  value  (r) 
comparisons  for  the  WSMR  site  10  m  stable  case  WOCSS  and  OB¬ 
WIND  (u„)  windfieid  outputs  where  n  is  the  number  of  surface  ob¬ 
servations  used  to  initialize  model  runs  with  a  250  m  grid. 


Component  Sources  Statistic  OBWIND  WOCSS 

U5-U5 

r 

1.000000 

1.000000 

U4-U5 

r 

0.107765 

1.000000 

U3-U5 

r 

0.121012 

1.000000 

uru5 

r 

0.099265 

1.000000 

Ul'Uj 

r 

0.280620 

1.000000 
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3.2  NTC  Results 

The  NTC  site  terrain  exhibited  greater  variance  than  any  of  the  other  3  sites 
incorporated  in  this  study  (see  Table  5  and  Fig.  10).  Therefore,  NTC  model  re¬ 
sults  represent  performance  characteristics  in  the  most  complex  terrain  available, 
the  opposite  extreme  of  the  WSMR  site  results. 

The  influence  of  change  in  terrain  complexity  is  reflected  in  comparing  some 
WSMR  and  NTC  results.  For  example,  under  stable  conditions  10%  (u  com¬ 
ponent)  and  78%  (v  component)  lower  average  correlation  coefficients  were  ob¬ 
served  in  comparing  stable  WSMR  and  NTC  results  (Tables  9-10  and  15-16). 
Conversely,  other  WSMR-NTC  result  comparisons  did  not  reflect  expected  re¬ 
sponses  to  terrain  complexity.  For  example,  under  stable  conditions  a  combined 
-17%  u  and  v  component  change  in  average  root  mean  square  differences  re¬ 
sulted  when  respective  WSMR  root  mean  sqaure  differences  were  compared  to 
NTC  root  mean  square  differences.  Also,  a  +52%  combined  (u  and  v  component 
average  agreement  index  value  change  between  the  WSMR  and  NTC  results  was 
unexpected.  This  example  is  more  reaffirmation  that  studies  should  not  rely  on 
single  statistical  indices  to  characterize  complex  associations. 

The  influence  of  change  in  terrain  complexity  was  more  consistently  reflected 
in  comparisons  between  the  WSMR  and  NTC  unstable  results.  For  example, 
under  unstable  conditions  Table  17  and  18  indicate  a  +44%  change  took  place  in 
the  combined  u  and  v  component  average  root  mean  square  differences  between 
respective  WSMR  and  NTC  results  (Tables  11-12  and  17-18);  lower  average  com¬ 
bined  agreement  index  values  (-75%),  as  well  as  -50%  lower  combined  average 
correlation  coefficients,  were  also  observed  in  comparing  WSMR  and  NTC  un¬ 
stable  case  results. 

Significant  differences  also  exist  between  the  NTC  stable  and  unstable  results. 
For  example,  average  combined  stable  case  agreement  index  values  were  29% 
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higher  than  the  unstable  case  agreement  index  values  (see  Tables  15-16  and 
17-18  for  respective  listings  of  u  and  v  component  agreement  index  values). 
Combined  u  and  v  component  average  stable  case  root  mean  square  differences 
and  correlation  coefficients  were  also,  respectively,  169%  and  326%  lower  than 
in  the  NTC  unstable  case,  which  suggests  that  it  is  important  to  consider  both 
stable  and  unstable  case  scenarios  in  high  resolution  simulation  applications. 

As  in  the  WSMR  results,  often  more  tangible  insight  into  study  results  was  ob¬ 
tained  from  analysis  of  graphical  spatial  depictions  of  model  outputs.  For  exam¬ 
ple,  Fig.  41  illustrates  the  basically  monotonic  OBWIND  model  wind  speed  sur¬ 
face  output  under  stable  conditions  at  25  m  AGL.  Figure  42  illustrates  WOCSS 
model  results  at  the  same  elevation  AGL  and  under  the  same  atmospheric  condi¬ 
tions.  Note  that  the  majority  of  the  WOCSS  25  m  level  solution  is  flow  restricted 
space,  which  is  a  dramatic  change  from  the  WOCSS  response  to  WSMR  terrain 
at  25  m  AGL  (see  Fig.  34).  Figure  43  illustrates  the  vector  differences  between 
the  25  m  AGL  stable  case  OBWIND  and  WOCSS  model  outputs  presented  in 
Figs.  41  and  42  respectively. 

Figure  44  demonstrates  that,  unlike  the  WSMR  unstable  case  results,  the 
WOCSS  model  produces  significant  unstable  case  flow  restricted  spaces  at  100 
m  AGL  using  NTC  site  terrain.  In  conjunction  with  Figs.  41-43,  Fig.  44  also 
demonstrates  the  WOCSS  model’s  dramatic  response  to  increases  in  terrain  com¬ 
plexity  (when  compared  to  the  WSMR  Fig.  34  results). 
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Table  15.  Summary  u  component  windfield  statistics  for  NTC  site  10, 
25,  and  100  m  stable  case  OBWIND  (u0)  and  WOCSS  (u^,)  windfield 
outputs  (rms=Root  Mean  Square  Difference,  agi=Agreement  Index, 
and  r=Correlation  Coefficient). 


Component  Sources  Statistic  100m  grid  250m  grid 

1000m  grid 

u0-u„, 

rms  10m 

0.565917 

0.568957 

0.581670 

u0-u„, 

rms  25m 

0.313170 

0.315G09 

0.353918 

Uq-Uu, 

rms  100m 

0.698934 

0.702900 

0.765826 

Uo-Uu, 

agi  10m 

0.406613 

0.384603 

0.281691 

u0-uw 

agi  25m 

0.323927 

0.301601 

0.221548 

Up-Uu, 

0.198802 

0.187169 

0.129889 

Up-Uu, 

r  10m 

0.011039 

0.011241 

-0.004228 

Up-uw 

r  25m 

0.139170 

0.160142 

0.206624 

Up-Uu, 

r  100m 

mmmum 

0.211968 

0.206624 

Table  16.  Summary  v  component  windfield  statistics  for  NTC  site  10, 
25,  and  100  m  stable  case  OBWIND  (v0)  and  WOCSS  (v„,)  windfield 
outputs  (rms=Root  Mean  Square  Difference,  agi=Agreement  Index, 
and  r=Correlation  Coefficient). 


Component  Sources  |  Statistic  }  100m  grid  (  250m  grid  |  1000m  grid 


rms  10m 


rms  25m 


rms  100m 


0.583754 


1.096843 


0.586630 


0.799898 


1.118825 


0.592768 


0.817913 


1.185041 


agi  10m 


25m 


0.075915 


0.028634 


0.049379 


0.018969 


0.044610 


0.015383 


0.018897 


0.116267 


-0.010592 


0.016418 


.095321 


-0.008645 


-0.183402 


0.025160 


-0.183583 


-0.093201 


r  10m 
r  25m 
r  100m 
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Table  17.  Summary  u  component  windfield  statistics  for  NTC  site  10, 
25,  and  100  m  unstable  case  OBW1ND  (u„)  and  WOCSS  (uw)  wind- 
field  outputs  (rms=Root  Mean  Square  Difference,  agi=Agreement 
Index,  and  r=Correlation  Coefficient). 


Component  Sources  Statistic  100m  grid  250m  grid 

1000m  grid 

Uo-Uu, 

rms  10m 

0.548071 

0.558534 

0.572248 

Ua-U,,, 

rms  25m 

0.681466 

0.702550 

0.733199 

u„-uw 

rms  100m 

2.114018 

2.755057 

u„-uw 

agi  10m 

0.324205 

0.296959 

0.198717 

u0-u„, 

0.147062 

0.144025 

0.157402 

Ua-Uu, 

0.061410 

-  -  ■  ■  -■ 

0.060150 

0.052337 

u0-uw 

r  10m 

0.009522 

-0.008320 

u0-uw 

r  25m 

0.201323 

0.237528 

0.331363 

Ua-Uw 

r  100m 

0.072375 

0.160770 

■riWirrm— 1 

Table  18.  Summary  v  component  windfield  statistics  for  NTC  site  10, 
25,  and  100  m  unstable  case  OBWIND  (v0)  and  WOCSS  (v„,)  wind- 
field  outputs  (rms=Root  Mean  Square  Difference,  agi=Agreement 
Index,  and  r=Correlation  Coefficient). 


Component  Sources  Statistic  |  100m  grid  250m  grid  ]  1000m  grid 


V0-Vw 

rms  10m 

1.507580 

1.515548 

1.514612 

v0-vw 

rms  25m 

mxmvm 

2.252696 

2.228796 

v0-vw 

rms  100m 

3.563609 

4.278866 

Vo-Vw 

agi  10m 

0.067367 

0.054279 

0.019853 

v„-vw 

agi  25m 

0.035105 

0.024261 

0.051423 

v0-vw 

FffWTiTi'ffW 

0.034248 

0.049633 

0.162873 

v.-v. 


v„-v„ 


r  10m 
r  25m 
r  100m 
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0.009106 

0.364159 

0.262189 


0.011506 

0.170301 

0.171201 


0.027654 

0.019495 

0.752704 
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grid  cells  (81  X  81 ).  NTC  site  terrain  elevation  is  plotted  below  the  speed  surface  plot  (north  is  located 
at  the  top  labeled  edge  of  the  plotted  regions).  Minimum,  average,  and  maximum  speed  values  are, 
respectively,  0.705, 0.754,  and  0.942  m/s,  where  speed  increases  along  the  vertical  plot  axis. 


8 


82 


Fiaure  44.  WOCSS 100  m  stable  case  wind  speed  plot  for  the  NTC  site  using  250  X  250  square  m 
arid  cells  (81  X  81).  NTC  site  terrain  elevation  is  plotted  below  the  speed  surface  plot  (north  is  located 
at  the  top  labeled  edge  of  the  plotted  regions).  Minimum,  average,  and  maximum  speed  values  are, 
respectively,  0.000, 1.141,  and  3.344  m/s,  where  speed  increases  along  the  vertical  plot  axis. 


3.3  Project  WIND  Results 

The  Project  WIND  site  terrain  exhibited  intermediate  variance  compared  to  the 
NTC  and  WSMR  sites  (see  Table  5  and  Fig.  10).  Therefore,  Project  WIND 
model  results  represent  performance  characteristics  of  intermediate  terrain  com¬ 
plexity,  and  the  influence  of  change  in  terrain  complexity  is  reflected  in  comparing 
some  WSMR  and  Project  WIND  results.  For  example,  under  stable  conditions 
a  -47%  combined  u  and  v  component  change  in  average  correlation  coefficients 
was  observed  in  comparing  stable  WSMR  and  Project  WIND  results  (Tables  9-10 
and  19-20). 

However,  as  in  the  WSMR-NTC  stable  case  comparisons,  other  results  did 
not  reflect  expected  responses  to  terrain  complexity.  For  example,  under  stable 
conditions  a  combined  -28%  u  and  v  component  change  in  average  root  mean 
square  differences  resulted  when  respective  WSMR  root  mean  square  differences 
were  compared  to  corresponding  Project  WIND  root  mean  square  differences. 
Also,  a  +81%  combined  (u  and  v)  component  average  agreement  index  value 
change  between  the  WSMR  and  Project  WIND  results  was  unexpected,  and  not 
well  correlated  with  changes  in  terrain. 

Mixed  results  were  also  obtained  in  WSMR-Project  WIND  unstable  case 
comparisons.  For  example,  under  unstable  conditions  Table  19  and  20  indicate^ 
that  average  combined  root  mean  square  differences  between  WSMR  (Tables  1 1- 
12)  and  Project  WIND  results  decreased  by  22%,  average  combined  agreement 
index  values  also  increased  by  62%,  and  average  combined  correlation  coefficients 
decreased  by  20%. 
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Noteworthy  differences  also  exist  between  Project  WIND  stable  and  unstable 
results.  For  example,  combined  average  Project  WIND  agreement  index  values 
were  50%  higher  in  unstable  case  (Tables  19-20  and  21-22).  Combined  average 
u  and  v  component  unstable  case  root  mean  square  differences  and  correlation 
coefficients  were  also  respectively  67%  and  169%  greater  than  in  the  Project 
WIND  stable  case. 

The  improved  agreement  index  and  correlation  coefficient  values  in  the  Project 
WIND  unstable  cases  are  attributed  largely  to  the  decrease  in  flow  restricted 
space  predicted  by  the  WOCSS  model  under  unstable  conditions  (e.g.,  in  the 
WSME  unstable  case  there  was  no  predicted  flow  restriction  at  10,  25,  or  100 
m  AGL  (see  Fig.  40)).  This  response  is  largely  due  to  the  increased  amount  of 
kinetic  energy  available  in  unstable  case  initialization  data  (see  Tables  A5-A8  in 
Appendix  A). 

As  in  the  WSMR  and  NTC  cases,  graphical  spatial  depictions  revealed  many 
other  interesting  results.  For  example,  Fig.  45  illustrates  a  10  m  AGL  stable 
case  OBWIND  speed  surface  output.  Note  the  speed  maximum  in  the  southeast 
corner  of  the  plot;  it  is  associated  with  higher  velocity  initialization  data  located 
in  the  vicinity  of  station  S-4  (see  Table  A7.1  in  Appendix  A).  Conversely,  Fig. 
46  illustrates  the  WOCSS  model  output  at  the  same  elevation  AGL  and  under 
the  same  conditions.  Most  of  the  WOCSS  solution  at  this  level  is  flow  restricted 
space,  except  for  several  regions  largely  concentrated  in  the  southwest  corner  of 
the  Project  WIND  study  area,  which  is  associated  with  a  low  level  plane  region. 

Figure  47  depicts  the  25  m  AGL  OBWIND  stable  case  output,  and  it  em¬ 
phasizes  a  speed  maximum  region  in  the  northeast  portion  of  the  Project  WIND 
study  area.  However,  the  corresponding  WOCSS  model  output  (Fig.  48)  presents 
a  very  different  prediction  solution:  the  OBWIND  speed  maximum  area  is  largely 
predicted  to  be  a  flow  restricted  space  by  the  WOCSS  model  (also  illustrated  in 
the  25  m  vector  difference  plot,  Fig,  49). 
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Table  19.  Summary  u  component  windfield  statistics  for  Project 
WIND  site  10, 25,  and  100  m  stable  case  OBWIND  (u0)  and  WOCSS 
(u,,,)  windfield  outputs  (rms=Root  Mean  Square  Difference,  agi= Agreement 
Index,  and  r=Correlation  Coefficient). 


Component  Sources  |  Statistic  j  100m  grid  j  250m  grid  |  1000m  grid 


rms  10m  I  0.646944 


rms  25m  0.443715 


rms  100m  0.531459 


0.646271 


.438451 


0.438451 


0.641437 


0.412346 


0.522843 


agi  10m 


agi  25m 


0.299979 


0.241593 


0.217753 


0.284134 


0.230301 


0.208817 


0.211191 


0.192204 


0.178661 


r  10m 


r  25m 


r  100m 


-0.059478 


-0.014800 


-0.082292 


0.048671 


0.038804 


Table  20.  Summary  v  component  windfield  statistics  for  Project 
WIND  site  10,  25,  and  100  m  stable  case  OBWIND  (v,)  and  WOCSS 
(vw)  windfield  outputs  (rms=Root  Mean  Square  Difference,  agi=Agreement 
Index,  and  r=Correlation  Coefficient). 


Component  Sources  |  Statistic  j  100m  grid  250m  grid  1000m  grid 


rms  10m 


rms  25m 


rms  100m 


.709004 


.500635 


0.708661 


0.499107 


0.844008 


0.705502 


0.490054 


0.853185 


agi  10m 


EM'ICTI 


0.291156 


0.202086 


0.151989 


0.279760 


0.222518 


0.179358 


.232102 


0.339326 


0.322698 


r  10m 


r  25m 


r  100m 


.114776 


0.219427 


0.152517 


0.113705 


0.269325 


.198306 


0.098359 


0.278671 


0.600944 
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Table  21.  Summary  u  component  windfield  statistics  for  Project 
WIND  site  10,  25,  and  100  m  unstable  case  OBWIND  (u0)  and 
WOCSS  (u*)  windfield  outputs  (rms=Root  Mean  Square  Difference, 
agi=Agreement  Index,  and  r=Corre!ation  Coefficient). 


Component  Sources  Statistic  100m  grid  250m  grid  1000m  grid 

Uo-U* 

rms  10m 

0.529500 

0.528440 

0.525441 

u0-uw 

rms  25m 

0.525855 

0.520805 

0.515802 

u0-uw 

rms  100m 

1.020713 

1.022311 

1.041712 

u0-u„, 

agi  10m 

0.399459 

0.408845 

0.437897 

Uo-Uu, 

agi  25m 

0.417741 

0.425086 

0.474907 

u0-uw 

agi  100m 

0.271813 

0.270955 

0.310996 

u0-uw 

r  10m 

0.210885 

0.090427 

u„*uw 

r  25m 

0.147367 

0.224680 

0.361363 

u0-uw 

r  100m 

0.097640 

0.186825 

0.444199 

Table  22.  Summary  v  component  windfield  statistics  for  Project 
WIND  site  10,  25,  and  100  m  unstable  case  OBWIND  (v0)  and 
WOCSS  (vw)  windfield  outputs  (rms=Root  Mean  Square  Difference, 
agi = Agreement  Index,  and  r=Correlation  Coefficient). 


Component  Sources  |  Statistic  100m  grid  |  250m  grid  |  1000m  grid  | 

v„-vw 

rms  10m 

0.864071 

0.877027 

0.928930 

v0*vw 

rms  25m 

0.890080 

0.876337 

0.763384 

v0-vw 

rms  100m 

2.160367 

2.182069 

2.204209 

v0-vw 

agi  10m 

0.395331 

0.358020  ' 

0.198166 

v0*vw 

iU-li&Ji 

0.380815 

0.425353 

0.630086 

v„-vw 

agi  100m 

0.134768 

0.169963 

0.334250 

v„-vw 

r  10m 

0.280300 

mmmim 

0.299044 

v0-vw 

0.132761 

0.488362 

Vo-Vw 

r  100m 

0.031702 

0.477423 

0.848978 
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3.4  Czechoslovakian  Results 


The  Czechoslovakian  site  terrain  had  intermediate  variance  compared  to  the  NTC 
and  WSMR  sites  (see  Table  5  and  Fig.  10).  Therefore,  the  Czechoslovakian 
model  results,  similiar  to  Project  WIND  Results,  also  represent  performance 
characteristics  in  transitional  terrain  complexity.  However,  the  Czechoslovakian 
results  represent  a  completely  different  set  of  initialization  conditions:  these  re¬ 
sults  were  produced  under  sparse  (single  point),  user  specified  initialization  con¬ 
ditions,  where  no  actual  observation  data  was  available.  This  exercise  was  pri¬ 
marily  implemented  to  study  sparse  initialization  behavior  of  the  WOCSS  model; 
therefore,  the  following  results  only  address  WOCSS  model  performance. 

An  extreme  condition  of  sparse  data  initialization,  typical  of  the  oversimplified 
conditions  presently  used  to  initialize  various  atmospheric  transport  processes  in 
a  number  of  U.S.  Army  applications,  was  employed  to  initialize  Czech.  WOCSS 
model  runs:  only  one  .vind  speed  and  direction  profile  was  used  to  initialize  each 
WOCSS  model  run.  The  same  wind  speed  and  direction  profile  was  used  to 
initialize  the  WOCSS  model  from  2  separate  locations  having  the  same  elevation 
(see  Appendix  A).  A  WOCSS  model  run  was  executed  for  each  separate  location 
(stations  S-l  and  S-2)  under  exactly  the  same  stable  and  unstable  conditions, 
which  produced  the  following  inner  WOCSS  model  comparison  results. 

For  example,  Tables  23-24  and  25-26  are  the  respective  u  and  v  component 
summary  statistics  for  the  stable  and  unstable  case  comparisons  of  stations  S-l 
and  S-2  (see  Appendix  A)  WOCSS  model  results.  Results  presented  in  Tables 
23-26  are  significantly  different  from  the  OBWIND-WOCSS  summary  statistics 
for  the  WSMR,  NTC,  and  Project  WIND  sites  (Tables  9-12, 15-18,  and  19-22,  re¬ 
spectively).  Yet  root  mean  square  differences  are  significantly  lower  for  the  Czech, 
inner  WOCSS  model  comparisons,  suggesting  little  difference  exists  between  S-l 
and  S-2  results;  this  result  is  complemented  by  relatively  high  agreement  index 
and  correlation  coefficient  values  (consistently  at  25  and  100  m  AGL). 
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However,  despite  the  significant  level  of  agreement  indicated  by  the  S-l/S-2 
WOCSS  results  in  Tables  23-26,  agreement  and  correlation  is  not  consistent  (e.g., 
see  Table  26)  and  suggests  that  important,  layer  selective  differences  can  result 
when  initializing  a  complex  terrain  wind  model  with  exactly  the  same  data  from 
different  points  over  the  simulated  terrain  space.  This  result  demonstrates  one 
of  the  flaws  in  assuming  a  single  mean  wind  speed  and  direction  over  a  complex 
terrain  space;  in  other  words,  assumptions  of  oversimplified  mean  wind  speed 
and  direction  over  a  complex  terrain  space  are  functionally  dependent  on  the 
initialization  location  within  the  terrain  space. 

Further  analysis  of  these  results  suggest  that  even  when  consideration  of  the 
functional  dependence  on  initialization  location  is  made,  the  assumption  of  using 
a  mean  wind  speed  and  direction  over  a  complex  terrain  space  is  unrealistic.  For 
example,  Figs.  50  and  51  represent  the  spatial  variation  in  10  and  25  m  AGL 
WOCSS  speed  surfaces  for  station  S-l  (see  Appendix  A)  under  stable  conditions. 
Notice  the  dramatic  variation  in  flow  restricted  space  between  the  10  and  25  m 
flow  surfaces  in  Figs.  50  and  51. 

Figure  52  illustrates  the  WOCSS  vector  differences  between  stations  S-l  and 
S-2  under  stable  conditions;  significant  vector  differences  are  observed  over  a 
widely  distributed  area,  not  just  at  initialization  locations.  Figures  53  and  54 
also  illustrate  the  significant  WOCSS  variation  produced  at  25  m  AGL  when 
station  S-2  was  initialized,  respectively,  with  stable  and  unstable  data,  and  the 
differences  between  station  S-2  and  station  S-l  at  25  m  AGL  under  unstable 
conditions  (see  Appendix  A)  are  plotted  in  Fig.  55. 

Figure  56  illustrates  the  magnitude  of  the  wind  speeds  associated  with  the 
vector  differences  indicated  in  Fig.  55.  Both  Figs.  55  and  56  suggest  that 
unstable  case  differences  were  very  dependent  on  placement  of  initialization  data 
in  the  modeled  Czechoslovakian  space  (see  Appendix  A  for  locations  of  S-l  and  S- 
2).  Because  terrain  influence  on  the  unstable  case  differences  is  minor  (compared 
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to  the  effect  of  initialization  data  placement),  there  appears  to  be  little  physical 
significance  in  the  unstable  case  results.  However,  the  distribution  of  vector 
differences  in  the  stable  case  comparison  (see  Fig.  52),  where  average  wind 
speeds  at  25  m  AGL  were  approximately  80%  less  than  average  unstable  case 
wind  speeds,  appears  to  be  much  more  influenced  by  terrain  effects.  For  example, 
stable  case  speed  differences,  as  illustrated  in  Fig.  57,  appear  to  be  closely 
associated  with  the  occurrence  of  rugged  terrain  areas,  which  required  further 
analysis. 

For  example,  terrain  slope,  the  rate  at  which  terrain  elevation,  z,  changes  in 
the  horizontal 

dz  =  (£)  * + (|)  * 

and  the  change  in  slope 

d'z = (B) 4x1 + (0)  ^ 

were  solved  for  over  the  selected  20  X  20  km2  Czechoslovakian  terrain  in  an  at¬ 
tempt  to  infer  physical  significance  in  the  ^table  case  result  differences  illustrated 
in  Figs.  52  and  57.  The  dz  and  dPz  values  were  translated,  respectively,  into  0 
and  d9  elevation  angles  above  the  x  -  y  tangent  plane. 

Figure  58(a)  and  58(b)  illustrate,  respectively,  the  0  and  d9  solutions  for  the 
Cech.  site,  which  agree  well  with  the  original  terrain  elevation  data  depicted  in 
Fig.  7(b)  -  for  example,  low  slope,  plain  areas  within  the  Cechoslovakian  terrain 
site  are  inversely  correlated  with  the  frequency  of  occurrence  of  plotted  contours 
in  Fig.  58. 


(15) 

(16) 
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In  analyzing  where  stable  wind  speed  differences  were  >  1.0  m/s,  it  was  dis¬ 
covered  that  100%  of  these  differences  occurred  between  the  475-525  in  terrain 
elevation  zones  within  the  20  x  20  km2  Czech,  site.  These  zones  were  repre¬ 
sentative  of  the  Czech,  site’s  mean  elevation  (e.g.,  the  niininimn,  mean,  and 
maximum  terrain  elevation  ASL  was,  respectively,  374.0,  188.0,  and  710.0  in). 
Analysis  was  then  focused  on  the  correlation  between  this  occurrence  of  wind 
speed  differences  >  1.0,  6,  and  d9  (see  Table  27). 

Table  28  indicates  that  approximately  100%  of  these  speed  differences,  >  1 .0 
m/s  between  the  475-525  m  ASL  terrain  zones,  occurred  along  relatively  low 
to  moderate  slopes  (0°  —  12°).  Table  29  presents  further  analysis  results  that 
indicated  a  significant  amount  (57%)  of  the  wind  speed  differences  >  1.0  m/s 
also  occurred  in  the  vicinity  of  moderate  changes  in  slope  (4°  -  8°). 

Figure  59  illustrates  the  widely  distributed  and  selective  correlation  between 
0,  d0 ,  and  wind  speed  differences  >  0.25  m/s.  As  'Fables  28  and  2!)  suggesi. 
significant  wind  speed  differences  do  not  generally  coincide  with  steepest  slopes 
or  changes  in  slope,  but  Fig.  59  does  show  that  significant  wind  differences  arc 
closely  coupled  to  areas  where  terrain  slope  varies  frequently.  Figures  57  and  5!). 
therefore,  appear  to  illustrate  a  physically  significant  Czech,  site  WOC-SS  model 
response  under  the  given  sparse  data  stable  conditions,  where  wind  speeds  were 
3.0  m/s;  and  this  particular  WOCSS  model  result  demonstrates  the  unrealistic 
value  in  assuming  a  single  wind  speed  and  direction  to  be  i -'present ative,  for 
example,  in  high  resolution  atmospheric  transport  applications  over  20  km  x  20 
km  moderately  complex  terrain  areas. 
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Table  23.  Summary  u  component  windfield  statistics  fwi  Cz«vU. 
site  10,  25,  and  100  m  stable  case  WOCSS  (u,lajt<mi)  and  WOCSS 
(u^ioiionj)  windfield  outputs  (rms=Root  Mean  Square  Difference,  agi 
=  Agreement  Index,  and  r=Correlation  Coefficient). 


Component  Sources  Statistic  250m  grid 

U,la.l*U4<a.2 

rms  10m 

0.002996 

rms  25m 

0.086433 

rms  100m 

0.029136 

UJl«.l'u»ia.2 

agi  10m 

0.999565 

U»ta.l’Ufta.2 

agi  25m 

0.959164 

tt»(a.l'Uf{0,3 

0.896098 

Uita.l*Uj*«.2 

r  10m 

0.999193 

Uj<a.|*UJ<a.2 

r  25m 

0.920526 

U*ia.l*u»<a.2 

r  100m 

0.809951 

Table  24.  Summary  v  component  windfield  statistics  for  Czech, 
site  10,  25,  and  100  m  stable  case  WOCSS  (v,fa<tonl)  and  WOCSS 
(v,<a«.ona)  windfield  outputs  (rms=Root  Mean  Square  Difference,  agi 
=  Agreement  Index,  and  r=Correlation  Coefficient). 


|  Component  Sources  Statistic  250m  grid  ( 

V,ta.l-V,ja.2 

rms  10m 

0.010131 

rms  25m 

0.222958 

rms  100m 

0.078130 

v*ta.l’v»f».2 

KS9XI2I 

0.999590 

agi  25m 

0.973557 

agi  100m 

0.835209 

V||j.rV|(a.3 

r  10m 

0.999179 

r  25m 

0.948924 

v»to.rv*i«.2 

r  100  m 

0.707207 
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Table  25.  Summary  u  component  windfield  statistics  for  Czech,  site 
10,  25,  and  100  m  unstable  case  WOCSS  (uJ<atloni)  and  WOCSS 
(Ujtaiionz)  windfield  output?  (rms=Root  Mean  Square  Difference,  agi 
=  Agreement  Index,  and  r=Correlation  Coefficient). 


|  Component  Sources  Statistic  250m  grid 

Uiii.rU,ii.i 

rms  10m 

0.472591 

U,|a.i-U*<a.2 

rms  25m 

0.047367 

U,ia.l‘U,|a.2 

rms  100m 

0.013240  1 

UJta.l*Uj|a.2 

agi  10m 

0.002275 

Uwa.l‘Uila.2 

agi  25m 

0.692139 

U»ta.l'Uj|a.2 

agi  100m 

0.976293 

Uita.l'll,«a.2 

r  10m 

-0.098538 

U»«a.l'Uita.2 

r  25m 

0.481701 

^ili.r^ila.2 

r  100m 

0.953453 

Table  26.  Summary  v  component  windfield  statistics  for  Czech,  site 
10,  25,  and  100  m  unstable  case  WOCSS  (v^iumi)  and  WOCSS 
(Vajationj)  windfield  outputs  (rms=Root  Mean  Square  Difference,  agi 
=  Agreement  Index,  and  r=Correlation  Coefficient). 


|  Component  Sources  Statistic  250m  grid  | 

V»i».rV|ia.2 

rms  10m 

1.611782 

V»|».rV|la.] 

rms  25m 

0.102231 

^»la.l‘V«ta,2 

rms  100m 

0.082092 

V*«a.l-Vj*a.2 

agi  10m 

^iii.r^iii.] 

agi  25m 

0.829984 

^*ta.l"V*ta.2 

agi  100m 

0.938916 

V>U.1-V  ill.] 


^ili.l*V|li.] 


r  10m 
r  25m 
r  100m 


-0.020724 

0.680182 

0.881847 
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Figure  50.  WOCSS  Station  1  stable  10  m  case  wind  speed  plot  for  the  Czech,  site  using  250  X  250 
square  m  grid  cells  (81  X  81).  Czech,  site  terrain  elevation  is  plotted  below  the  speed  surface  plot  (north 
is  located  at  the  top  labeled  edge  of  the  plotted  regions).  Minimum,  average,  and  maximum  speed  values 
are,  respectively,  0.000, 0.845,  and  2.202  m/s,  where  speed  increases  along  the  vertical  axis. 
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Figure  51 .  WOCSS  Station  1  stable  25  m  case  wind  speed  plot  for  the  Czech,  site  using  250  X  250 
square  m  grid  ceils  (81  X  81 ).  Czech,  site  terrain  elevation  is  plotted  below  the  speed  surface  plot  (north 
is  located  at  the  top  labeled  edge  of  the  plotted  regions).  Minimum,  average,  and  maximum  speed  values 
are,  respectively,  0.000, 0.812,  and  2.202  m/s,  where  speed  increases  along  the  vertical  axis. 


WOCSS  STA1 -STA2  VECTOR  DIFFERENCE  25M  STABLE  25QM  GRIP 


Figure  52.  WOCSS  Station  1-  Station  2  stable  25  m  case  wind  vector  difference  plot  for  the  Czech, 
site  using  250  X  250  square  m  grid  cells  (81  X  81 ).  Czech,  site  terrain  elevation  is  plotted  below  the 
speed  surface  plot  (north  is  located  at  the  top  labeled  edge  of  the  plotted  regions).  Magnitude  of 
average  vector  differences  -  0.066  m/s,  and  magnitude  of  maximum  vector  difference  -1.811  m/s. 
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Figure  53.  WOCSS  Station  2  stable  25  m  case  wind  speed  plot  for  the  Czech,  site  using  250  X  250 
square  m  grid  cells  (81  X  81).  Czech,  site  terrain  elevation  is  plotted  below  the  speed  surface  plot  (north 
is  located  at  the  top  labeled  edge  of  the  plotted  regions).  Minimum,  average,  and  maximum  speed  values 
are,  respectively,  0.000, 0.812,  and  2.196  m/s,  where  speed  increases  along  the  vertical  axis. 
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Figure  54.  WOCSS  Station  2  unstable  25  m  case  wind  speed  plot  for  the  Czech,  site  using  250  X  250 
square  m  grid  cells  (81  X  81).  Czech,  site  terrain  elevation  is  plotted  below  the  speed  surface  plot  (north 
is  located  at  the  top  labeled  edge  of  the  plotted  regions).  Minimum,  average,  and  maximum  speed 
values  are,  respectively,  4.153, 4.191 ,  and  5.533  m/s,  where  speed  increases  along  the  vertical  axis. 


WOCSS  STA1-STA2  VECTOR  DIFFERENCE  25M  UNSTABLE  25QM  GRID 
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Figure  56.  WOCSS  Station  1-  Station  2  unstable  25  m  case  wind  sped  difference  plot  for  the  Czech, 
site  using  250  X  250  square  m  grid  cells  (81  X  81 ).  Czech,  site  terrain  elevation  is  plotted  below  the 
speed  surface  plot  (north  is  located  at  the  top  labeled  edge  of  the  plotted  regions).  Magnitude  of  average 
speed  differences  *  0.063  m/s,  and  magnitude  of  maximum  speed  difference  *  1 .292  m/s,  with  speed 
increasing  upward  along  the  vertical  plot  axis.. 
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Table  27.  Local  maximum  slope,  0,  and  changes  in  slope,  d$,  in 
degrees,  for  the  Czechoslovakian  20  km  x  20  kma  terrain  site  using 
100  m  grid  spacing. 


9 

uu 

MINIMUM 

0.008 

hkium 

MEAN 

KX£]j|| 

MAXIMUM 

17.0°  || 

Table  28.  Correlation  (in  percent)  of  9  and  the  occurrence  of  25  m 
AGL  stable  case  wind  speed  differences  greater  than  or  equal  to  1.0 
m/s  for  the  Czechoslovakian  20  km  x  20  km3  terrain  site. 


|  9  (in  degrees) 

wind  speed  differences  >  1.0  m/s 

F  49  % 

1  51  % 

umpfrareai 

'■■■HHilj 

Table  29.  Correlation  (in  percent)  of  d9  and  the  occurrence  of  25  m 
AGL  stable  case  wind  speed  differences  greater  than  or  equal  to  1.0 
m/s  for  the  Czechoslovakian  20  km  x  20  km3  terrain  site. 


wind  speed  differences  >1.0  m/s  | 

43  %  ”“1 

II  4®  <  d0  <  8° 

57  %  1 

f-  ~TK 
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4  CONCLUSIONS 


Under  the  light  wind  speeds  used  to  initialize  all  stable  and  unstable  case  runs 
(see  Appendix  A),  WSMR,  NTC,  and  Project  Wind  WOCSS  10  and  25  m  AGL 
results  generally  contributed  significant  supplementary  physical  insight  into  flow 
dynamics  that  were  missed  by  the  OB  WIND  model.  The  exaggerated  extent  of 
WOCSS  flow  restricted  space  at  10  m  AGL  for  Project  Wind  stable  conditions 
(Fig.  46),  and  for  NTC  25  m  AGL  stable  conditions  (Fig.  42),  suggest  that  the 
WOCSS  model  is  oversensitive  to  terrain  complexity  greater  than  that  at  the 
WSMR  site  (e.g.,  Figs.  25  and  34  demonstrate  a  more  reasonable,  moderated 
application  of  flow  restricted  space  prediction). 

However,  the  validity  of  WOCSS  model  responsiveness  to  variations  in  inter¬ 
mediate  to  very  complex  terrain  appears  to  improve  as  wind  speeds  are  increa'ed 
more  toward  neutral  conditions  (see  Table  6  and  Table  A9).  Under  the  light  wind 
speeds  considered  in  this  study  the  WOCSS  model  appears  to  predict  reasonably 
valid  flow  restricted  space  predictions  for  use  in  terrain  complexity  <  to  the 
WSMR  site  terrain  (see  Table  5  and  Fig.  10). 

Conversely,  the  WOCSS  model  has  some  serious  limitations  that  aie  never 
compensated  for.  For  example,  localized  horizontal  wind  direction  and  speed 
changes  (resulting  from  mechanical  forcing  that  occurs  in  nature  where  flow  sur¬ 
faces  intersect  with  terrain  obstacles)  are  not  calculated  (e.g.,  see  Fig.  24).  This 
limitation  severely  restricts  the  utility  of  unmodified  WOCSS  windfields  in  high 
resolution  atmospheric  transport  application  work. 

But  results  from  Table  6  and  Figs.  19-21  indicate  that  the  WOCSS  model 
demonstrates  significant  and  reasonable  responsiveness  to  typical  stable  and  un¬ 
stable  variations  in  boundary  layer  conditions,  as  they  are  observed  to  behave 
in  nature;  this  WOCSS  windfield  responsiveness  in  predicting  reasonable  strati¬ 
fied  variation  in  flow  restricted  spaces  could  be  incorporated  as  a  supplemental 
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initialization  improvement  for  use  in  another  model.  For  example,  initialization 
data  modified  to  incorporate  the  effects  of  flow  restricted  spaces,  predicted  by 
the  WOCSS  model,  could  then  be  used  to  initialize  a  basic  mass  consistency 
model,  such  as  OBWIND,  thereby  improving  the  overall  validity  of  composite 
model  outputs. 

Results  from  this  study,  therefore,  suggest  that  the  WOCSS  model  could  be 
selectively  employed  to  produce  improved  and  reasonable  estimations  of  terrain- 
windfield  coupling  effects  for  applications  that  have  not  previously  considered 
1)  the  dynamics  in  forcing  flow  over  mechanical  obstacles;  and  2)  the  influence 
of  boundary  layer  kinetic  energy  and  stability  effects  on  flow  fields,  flu'  results 
listed  in  Appendix  C  (Tables  C4  and  C5)  also  suggest  that  WOCSS  solutions 
should  not  present  significant  run-time  burdens  to  users  of  the  WOCSS  model. 
The  monotonic,  performance  of  WOCSS  interpolation  operations,  also  evident  in 
the  results  of  this  study,  are  an  algorithm  design  problem  that  could  be  restruc¬ 
tured  and  improved  in  the  future;  therefore,  this  performance  characteristic  is 
not  a  permanent  detraction  from  the  WOCSS  model’s  potential  validity. 

The  research  work  presented  here,  therefore,  demonstrates  that  the  WOCSS 
model  can  be  employed  to  simulate  complex  terrain  flow  dynamic  effects  better 
than  basic  mass  balancing/interpolatiou  schemes,  such  as  tin'  OBWIND  model. 
Measures  of  the  WOCSS  model's  performance  presented  in  this  report  suggest 
that,  although  it  does  not  simulate  dynamic  atmospheric  How  struct  lire  in  com¬ 
plex  terrain  with  exceptional  fidelity,  the  WOCSS  model  has  the  potential  to 
provide  improved  dynamic  flow  structure  approximations  for  a  variety  of  appli¬ 
cation  purposes  that  currently  incorporate  only  kinematic  solutions  or  that  use 
only  a  single  selected  mean  wind  speed  and  direction  over  a  complex  terrain  sur¬ 
face;  in  this  context,  the  WOCSS  model  could  provide  solutions  of  atmospheric 
flow  structures  that  have  a  significantly  improved  degree  of  validity. 


Potential  application  of  this  WOCSS  model  utility  could  be  incorporated,  for 
example,  at  the  National  Training  Center  (NTC)  in  California  where  the  U  S. 
Army  carries  out  low  level  helicopter  training  exercises  year-round.  Terrain  is 
very  rugged  at  the  NTC  (see  Figs.  5  and  10),  and  flight  forecast  support  is 
limited  there,  partially  because  the  NTC  occupies  a  relatively  data  sparse  area, 
and  the  availability  of  experienced  forecasters  permanently  stationed  at  the  NTC 
is  severely  limited  thus  making  it  an  ideal  location  for  reliable,  complex  terrain 
wind  model  applications.  WOCSS  model  guidance  could  be  used,  by  local  NTC 
forecasters  to  supplement  their  present  low  level  flight  hazard  analysis  schemes, 
thereby  improving  accuracy  of  hazard  area  predictions,  which  would  improve 
upon  the  safety  of  flight  operations  at  the  NTC. 

Other  priority  applications  that  could  benefit  from  incorporating  more  valid 
wind  field  initialization  include  the  U.S.  Army's  Combined  Obscuration  Model  for 
Battlefield- Induced  Contaminants  (COMBIC),  developed  by  Hoock  et  al.  (1987). 
The  effectiveness  of  applications  such  as  COMBIC  are  very  dependent  on  the  reli¬ 
ability  of  input  environmental  factors,  such  as  wind  direction  and  speed.  Selective 
use  of  the  WOCSS  model  to  produce  inputs  for  complex  terrain  COMBIC  appli¬ 
cations  has  the  potential  to  significantly  improve  current  COMBIC  methodology, 
which  could  ultimately  benefit  other  models  that  employ  COMBIC  inputs. 

Although  difficult  to  quantify,  improved  decision  making  and  resource  effec¬ 
tiveness  could  be  realized  from  incorporating  the  WOCSS  model  into  application 
environments,  such  as  the  ones  referred  to  above.  These  improvements  will  re¬ 
sult  from  more  well  founded  decisions  made  by  users  of  decision  aid  applications 
that  require  valid  atmospheric  flow  structure  information,  even  under  data  sparse 
circumstances. 

Technical  results  from  this  study  suggest  that  250  m  x  250  m  grid  spacing 
is  a  practical  limit  in  reporting  on  the  high  resolution  performance  of  modeling 
and  simulation  results  in  20  km  x  20  km  spaces.  The  250- m  grid  spacing  is  also 


probably  acceptable  for  high  resolution  field  application  use  in  rural  settings, 
dependent  on  detailed  analysis  of  terrain  complexity  (Fig.  10).  Even  the  1000-m 
grid  spacing  could  be  used  if  execution  run-time,  and  available  storage  device 
bad  to  be  reduced,  but  there  would  be  significant  degradations  in  resolving  local 
flow  feature  detail  (e.g.,  compare  Figs.  31  and  32).  However,  grid  spacing  < 
100  m  is  recommended  for  urban  scale  applications  that  address  the  influence  of 
mechanical  forcing  induced  by  significant  urban  structures  (e.g.,  skyscraper  and 
industrial/commercial  complexes). 

Other  results  of  this  study  suggest  that  flow  layer  output  that  is  used  as  input 
to  other  applications  must  be  carefully  selected  because  of  the  WOCSS  model’s 
stratified  responsivenese  to  initial  conditions  (see  Figs.  19-21  and  Table  6).  For 
example,  the  selection  of  the  flow  surface  elevation  must  be  significantly  related 
to  mean  obscurant  transport  levels  in  COMBIC  applications,  which  are  largely 
buoyancy  considerations. 

Results  from  this  study  also  suggest  that  no  single  statistical  measure  used 
here  (e.g.,  root  mean  square  difference,  agreement  index,  or  correlation  coeffi¬ 
cients)  can  be  relied  upon  to  infer  a  complete  characterization  of  model  valid¬ 
ity,  especially  without  the  assistance  of  supplemental  graphic,  spatial  analyses. 
For  example,  in  this  study,  OBWIND-WOCSS  correlation  coefficients  generally 
tended  to  increase  as  grid  spacing  resolution  and  predicted  flow  restricted  space 
decreased  (e.g.,  as  distance  from  the  terrain  surface  increased  in  the  vertical). 
However,  agreement  index  values  produced  smaller  and  less  consistent  changes 
when  grid  space  resolution  decreased;  agreement  index  changes  were  generally 
more  responsive  to  predictions  of  flow  restricted  space  in  the  vertical  stratifi¬ 
cation  comparisons,  but  this  trend  was  also  inconsistent  in  most  of  the  study 
results. 
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APPENDIX  A 


Meteorological  Initialization  Data 


Surface  meteorological  observation  towers  that  provided  observation  data  in 
this  research  work  had  anemometer  wind  sensors  elevated  at  approximately  ten 
meters  above  ground  level.  Sensor  locations  at  the  WSMR,  NTC,  and  Project 
WIND  sites  are  indicated  in  Figs.  Al-3,  and  Tables  A 1-3.  Radiosonde  upper- 
air  observation  data  was  also  available  at  the  WSMR  and  Project  WIND  sites. 
Figure  A4  and  Table  A4  note  the  initialization  points  for  the  Czechoslovakian 
site,  for  which  no  actual  meteorological  data  was  available.  Tables  A5-8  list  the 
actual  surface  initialization  data  used  in  the  case  studies. 

Case  study  initialization  data  were  broken  down  into  stable  and  unstable  con¬ 
ditions.  Several  stability  criteria  were  used  to  analyze  all  observation  data  sets. 
For  example,  Pasquill’s  Stability  Categories  1  were  relied  upon  to  initially  classify 
the  stability  rating  of  each  surface  observation  set,  primarily  using  wind  speed 
and  time  of  day  (insolation  inferencing).  Table  A9  lists  the  Pasquill  Stability 
Category  scheme  that  was  used  to  carry  out  an  initial  classification.  Local  diur¬ 
nal  trends  in  pressure,  and  temperature  were  taken  into  consideration  to  further 
rate  each  observation  set’s  stability  ratings. 

Temperature  lapse  rate  in  the  vertical  was  also  considered  where  radiosonde 
data  was  available  (see  Figs  A5-A8).  Potential  temperature,  0,  is  plotted  as  a 
function  of  height  above  ground  level  in  Figs.  A5-A8  and  was  solved  for  using 
the  following  equation  (from  Hess,  1979): 


9  =  77— P 

Kmo’ 


(A  -  1) 


where  T  is  the  temperature  of  a  parcel  of  air  that  has  been  moved  adiabatically 
lArya,  1988. 
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to  a  new  pressure,  p.  The  k  term3  is  equivalent  to  R/cp,  where  R  represents  the 
specific  gas  constant,  and  Cp  is  the  specific  heat  capacity  at  constant  pressure. 

The  potential  temperature  lapse  rate,  gj,  in  Figs.  A5-A8  illustrate  the  vertical 
stability  classification  according  to  the  criteria  in  Table  A 10.  Actual  vertical 
temperature  profile  data  was  used  for  the  WSMR  and  Project  WIND  model 
runs.  These  data  were  then  modified  to  develop  representative  initialization  data 
for  the  NTC  and  Czechoslovakian  sites.  It  should  be  noted  that  the  increase  in 
non-linear  behavior  of  the  WSMR  vertical  temperature  profile,  in  Fig.  A5,  over 
that  of  the  Project  WIND  data  (Fig.  A7)  resulted  from  a  higher  frequency  of 
observation  data  being  collected  during  the  WSMR  radiosonde  ascents. 

The  actual  WSMR  and  Project  WIND  are  noted  in  Table  All.  These 
data  were  used  to  develop  vertical  lapse  rate  initialization  data  for  the  NTC  and 
Czechoslovakian  model  runs.  For  example,  the  observed  lapse  rates  for  WSMR 
and  Project  WIND  were  averaged  and  then  reduced  by  a  factor  of  50%.  The 
50%  reduction  was  implemented  to  reduce  the  generated  lapse  rates  to  more 
neutral  conditions,  a  major  assumption  in  the  application  of  a  logarithmic  veloc¬ 
ity  profile  equation,  which  is  discussed  after  the  temperature  profile  construction 
methodology  has  been  reviewed. 

Several  steps  were  used  to  construct  theoretical  vertical  temperature  profiles 
for  the  NTC  and  Czechoslovakian  sites.  First,  the  pressure  at  sea  level  was 
estimated  for  each  site  under  both  the  stable  and  unstable  thermal  conditions 
(noted  in  Tabes  A6  and  A8,  respectively).  This  required  extrapolating  site  tem¬ 
peratures  to  sea  level  using  the  standard  atmospheric  lapse  rate  of  6.5°C/km 
(Huschke,  1989). 

**  ha*  a  value  of  0.286  for  dry  air  (Hess,  1979). 
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Then  an  equation  of  state  was  used  to  calculate  a  representative  sea  level 
pressure  (Hess,  1979): 

p  =  pRT  (A  -  2) 

where  p  is  pressure,  p  is  density  (set  at  0.0012250  gm /cm3,  from  Huschke,  1979), 
R  is  a  specific  gas  constant  (2870.4  mb  cm3/gm  K  for  dry  air,  from  Riegel,  1974), 
and  T  is  the  temperture  (in  Kelvin). 

The  hypsometric  formula  was  then  used  to  solve  for  the  pressure  of  each 
observation  site  under  each  stability  condition  (from  Riegel,  1974): 

Z3-Zl  =  /Wr*^  (A -3) 

9  jp%  P 

where  zj  —  z\  is  the  difference  in  vertical  elevation  between  the  points  of  interest, 
Rd  is  the  specific  gas  constant  for  dry  air,  g  is  acceleration  due  to  gravity,  pi 
and  pi  are  the  pressures  at  points  Z\  and  z 2,  respectively,  and  T’  is  the  virtual 
temperature. 

T*  can  be  found  using  the  following  equation,  also  from  Riegal  (1974): 

T-  =  T  (1  +0.608 - ~~ — r- )  (A -4) 

\  p_(l_e)e/ 

where  T  is  the  observed  temperature,  e  is  the  ratio  of  the  mass  of  water  vapor 
to  the  mass  of  dry  air,  e  is  the  partial  pressure  due  to  water  vapor,  and  p  is  the 
total  pressure.  In  this  research  work,  only  dry  air  cases  are  considered;  therefore, 
both  e  and  e  are  zero.  As  a  result,  T*  «  T  is  a  valid  assumption,  and  T  can  be 

used  in  equation  A-3.  After  solving  equation  A-3  for  pj,  then  equation  A-l  can 

be  used  to  find  0  for  each  case  at  each  site.  From  this  point,  Table  All  can  be 
used  to  quickly  construct  vertical  temperature  profiles  (see  Figs.  A6  and  A8). 

Since  radiosonde  data  were  not  available  for  the  NTC  and  Czechoslovakian 
sites,  hypothetical  wind  initialization  data  were  developed  using  a  logarithmic 
velocity  profile  approximation  (Arya,  1988)  assuming  a  neutral  boundary  layer 

U  =  u]-\n-  (A -5) 
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where  U  is  the  mean  horizontal  wind  speed  as  a  function  of  friction  velocity, 
u*,  von  Karman’s  constant,  k ,  the  roughness  parameter,  20,  and  height  above  a 
surface,  z.  The  friction  velocity  term  is  defined  as  u*  =  ypf,  where  r0  is  surface 
shear  stress,  and  p  is  the  fluid  density. 

Variability  of  the  NTC  terrain  elevation  data  was  greater  t  han  any  other 
site  (see  Table  5).  Therefore,  an  approximation  for  z0  at  the  NTC  was  made 
(from  Arya,  1988),  and  this  value  was  scaled  to  the  fractional  variability  of  the 
Czechoslovakian  site’s  elevation  data  (see  Table  A 12).  The  u"  term  used  for  the 
NTC  vertical  wind  profde  was  solved  for  using  actual  observation  data  at  2=10 
m,  and  the  selected  value  for  z0.  Then  the  logarithmic  velocity  profile  equation 
was  employed  to  develop  the  vertical  wind  profile  initialization  data  for  the  N  I  C 
model  runs.  A  similar  approach  was  used  to  generate  initialization  profile  data 
for  the  Czechoslovakian  model  runs,  except  that  the  su.face  winds  were  user 
specified  instead  of  being  actual  observation  data  (see  Tables  A 12  and  A 13). 

To  develop  accurate  vertical  wind  direction  profiles  the  theoretical  Ekman 
spiral  was  employed  through  the  use  of  the  following  equations,  also  from  Hess 
(1979), 

ti  =  ug[l  -  J-y/MK")  cos(v///2A'  •  z))  (A  -  G) 

v  =  sin(/f/2A'  •  z)  (A  -  7) 

where  u  and  v  arc  wind  speed  components  in  the  x  and  y  directions  respectively, 
/  is  the  Coriolis  parameter  («  10-4s_1),  K  originated  as  a  eddy-exchange  coef¬ 
ficient  constant  (empirically  set  at  «  5  x  104 cm2s-1),  and  2  represents  vertical 
elevation  above  ground  level.  An  Ekman  boundary  layer  balances  coriolis,  pres¬ 
sure  gradient,  and  frictional  forces  with  height,  up  to  the  lowest,  atmospheric 
level  free  from  surface  effect  influences;  winds  at  and  above  this  level  are  referred 
to  as  geostrophic  and  are  represented  above  as  the  ug  and  vg  components  (Hess, 
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1979).  Plots  of  the  modeled  wind  profiles  used  in  the  NTC  and  Czecoslovakian 
simulation  runs  are  illustrated  in  Figs.  A6  and  A8,  respectively. 

In  summary,  the  averaged  observation  data  at  the  WSMR  and  Project  WIND 
site  had  slightly  stronger  wind  velocities  («  2.46m/s  at  10  m  above  ground  level) 
for  the  unstable  cases  versus  stable  cases.  The  unstable  cases  also  had  greater 
average  magnitudes  in  potential  temperature  lapse  rates  (s»  0.83°C)  within  the 
first  100  m  above  ground  level  (see  Figs.  A5  and  A7).  These  properties  of  the 
observed  stable  and  unstable  conditions  were  carried  over  in  the  construction 
of  missing  data  for  the  NTC  and  Czechoslovakian  sites,  but  in  these  cases  the 
properties  were  adjusted  toward  neutral  conditions  -  to  accommodate  the  use  of 
equations  A-5  -  A-7  (see  Figs  A6  and  A8). 

Part  of  the  meteorological  initialization  data  presented  here  was  acquired  from 
actual  field  observations.  The  WSMR  and  Project  WIND  sites  had  complete  sets 
of  10-meter  surface  observation  data  that  were  coupled  to  radiosonde  vertical 
profile  data.  The  NTC  site  illustrated  a  transition  toward  more  data  sparse 
field  conditions,  where  only  a  limited  number  of  10  meter  surface  observations 
were  available.  However,  using  classical  methods,  correlated  vertical  data  were 
developed  to  establish  complete  sets  of  reasonable  model  initialization  data. 

Then  transition  toward  the  worst  case  scenario,  in  terms  of  data  sparsity,  was 
established  for  the  Czechoslovakian  site;  no  observation  data  was  available  for 
initializing  model  runs  for  this  site.  However,  by  tailoring  properties  of  obser¬ 
vation  data  from  the  other  sites,  in  conjunction  with  the  selective  application 
of  classical  meteorological  theory,  reasonable  model  initialization  data  was  also 
developed  for  the  Czechoslovakian  site.  The  construction  of  this  hypothetical, 
user  specified,  initialization  data  illustrates  methods  that  could  be  used  to  im¬ 
prove  the  credibility  and  resolution  of  environmental  conditions  in  remote  site 
simulation  activities. 
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Figure  A1.  WSMR  surface  meteorological  observation  tower  locations  (S-1 
through  S-5),  and  the  radiosonde  observation  launch  site  (UA). 
UTM  coordinate  references  are  included. 


Table  Al.  WSMR  meteorological  observation  sensor  locations. 
Observation  Site  ||  UTM  Easting  ||  UTM  Northing  ||  Surface  Elevation  (m) 


392426.7 

390597.0 

392637.9 

391510.3 

391632.1 

391483.2 


3595988.4 

3588896.5 
3588404.9 
3585415.8 
3585256.1 

3589972.6 


1250.0 

1254.1 
1288.9 
1282.3 

1322.2 
1253.0 


12 


3932000 


Figure  A2.  NTC  surface  meteorological  observation  tower  locations  (S-1 

and  S-2).  No  radiosonde  sites  were  available  in  this  study  area 
UTM  coordinate  references  are  included. 


Table  A2.  NTC  meteorological  observation  sensor  locations. 


Observation  Site 

UTM  Easting  UTM  Northing  Surface  Elevation  (m) 

S-1 

539500.0 

3920200.0  1180.0 

S-2 

547800.0  ) 

3927500.0  1080.0 

ii5 


Figure  A3.  Project  WIND  surface  meteorological  observation  tower  locations 
(S-1  through  S-5).  The  radiosonde  observation  launch  site  was  co¬ 
located  at  tower  S-4's  location.  UTM  coordinate  references  given. 


Table  A3.  Project  WIND  meteorological  observation  sensor  locations. 


Observation  Site 

UTM  Easting 

UTM  Northing 

Surface  Elevation  (m) 

S-1 

561170.0 

4459350.0 

216.0 

S-2 

564400.0 

4455350.0 

137.0 

S-3 

577760.0 

4456650.0 

543.0 

S-4 

574500.0 

4451100.0 

253.0 

S-5 

577380.0 

4452180.0 

271.0  ^ 
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Figure  A4.  Czechoslovakian  surface  meteorological  initialization  points  (S-1 
and  S*2).  No  actual  surface  or  upper  air  data  was  used  for  this  study  area; 
meteorological  data  was  strictly  user  defined.  UTM  references  given. 


Table  A4.  Czechoslovakian  meteorological  data  initialization  points. 


Initialization  Points 

UTM  Easting  ||  UTM  Northing  |{  Surface  Elevation  (m) 

S-1 

342050.0  1  5503770.0  ||  450.0 

S-2 

342050.0  1  5488770.0  |  450.0 

127 


Table  A5.1  WSMR  surface  wind  observation  data  (15  minute  aver¬ 
ages)  under  diurnally  stable,  early  morning  atmospheric  conditions 
at  03:00  L  on  21  July  1991. 


Wind  Direction  (deg) 

Temp.  (C) 

049 

21.64 

041 

22.01 

051 

22.56 

056 

22.54 

080 

22.53 

055 

22.26 

Table  A5.2  WSMR  surface  wind  observation  data  during  unstable, 
late-afternoon  conditions  at  16:00  L  on  27  July  1991. 


Station  |j  Wind  Speed  (m/s)  |{  Wind  Direction  (deg)  ||  Temp.  (C) 


30.17 


30.18 


29.54 


30.79 


29.99 


30.13 


Table  A6.1  NTC  surface  wind  observation  data  (15  minute  averages) 
under  stable,  early  morning  conditions  at  01:00  L  on  24  February 
1988. 


Wind  Direction  (deg) 

Temp.  (C) 

173 

12.78 

120 

14.44 

147 

13.61 

Table  A6.2  NTC  surface  wind  observation  data  during  unstable,  mid¬ 
afternoon  conditions  at  15:00  L  on  22  February  1988. 


Station  U  Wind  Speed  (m/s)  I  Wind  Direction 


20.00 


20.00 


20.00 
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Table  A7.1  Project  WIND  surface  wind  observation  data  under  di- 
urnally  stable,  late  evening  atmospheric  conditions  at  23:00  L  on  25 
June  1985. 


irection 


s-1  1 

1.96 

345 

27  A0 

2J2 

341 

28.00 

■■ 

0.81 

300 

*** 

S-4 

2.64 

307 

1.86 

306 

28.00 

AVG 

1.88 

320 

27.84 

Table  A7.2  Project  WIND  surface  wind  observation  data  during  un¬ 
stable,  mid-afternoon  atmospheric  conditions  at  15:00  L  on  25  June 
1985. 


Station 

Wind  Speed  (m/s) 

Wind  Direction  (deg) 

Temp.  (C) 

S-1 

3.16 

169 

33.80 

S-2 

2.59 

152 

34.20 

5.02 

174 

31.70 

mm 

5.30 

212 

33.50 

iKsai 

5.66 

197 

33.80 

AVG 

4.35 

181 

33.40 

Table  A8.1  Hypothetical  surface  wind  initialization  data  for  the  Czechoslo¬ 
vakian  terrain  representative  of  stable  conditions. 


Station 

Wind  Speed  (m/s) 

Wind  Direction  (deg) 

Temp.  (C) 

S-1 

1.02 

315 

iHmi 

IIK99I 

1.02 

315 

Table  A8.2  Hypothetical  surface  wind  initialization  data  for  the  Czechoslo¬ 
vakian  terrain  representative  of  unstable  conditions. 


d  Direction  (deg)  Temp. 


20.00 


20.00 


Table  A9.  Pasquill’s  Stability  Categorization  scheme  3  used  to  ini¬ 
tially  classify  surface  meteorological  data  into  stable  and  unstable 
sets.  £xf=extreme,  Mod=moderate,  and  5/<=slight. 


Wind  Speed  (m/s) 

Mod-Slg  Daytime  Insolation 

Nighttime 

i  <2 

Ext-Mod  Unstable 

i  2 

Ext-Mod  Unstable 

Sit-Mod  Stable 

4 

Mod-Sit  Unstable 

Near  Neutral 

6 

Sit  Unstable-Near  Neutral 

Near  Neutral 

>6 

Sit  Unstable-Near  Neutral 

Near  Neutral 

Table  A10.  Potential  temperature  lapse  rate,  vertical  stability 
classification  scheme  4 


Lapse  Rate 

Stability  Rating 

dd/dz  <  0 

UNSTABLE 

dQJdz  =  0 

NEUTRAL 

dO/dz  >  0 

STABLE 

Table  All.  Potential  temperature  lapse  rate,  in  actual  observation 
data  selected  for  the  WSMR  and  Project  WIND  model  runs.  Here  dz 
represents  the  distance  between  ground  level  ( sfc )  and  100  m  above 
ground  level,  and  dO  is  the  difference  between  0ioo  m  and  0,jc. 


Site 

Stable  dOjdz 

Unstable  dOjdz 

WSMR 

+  1.85° 

-2.60° 

Project  WIND 

+2.95° 

-3.85° 

AVERAGE 

+2.40° 

-3.23° 

I 

I 


3Arya,  1988. 

4lless,  1979. 
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Elevation  Above  Ground  Level  (m)  Elevation  Above  Ground  Level  (m) 


Figure  A5.  WSMR  stable  (a)  and  unstable  (b)  meteorological  upper  air  initialization 

data  correlated  with  the  surface  observation  data  in  Tables  A5.1  and  A5.2, 
respectively. 
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Elevation  Above  Ground  Level  (m)  Elevation  Above  Ground  Level  (m) 


Figure  A6.  NTC  stable  (a)  and  unstable  (b)  modeled  upper  air  initialization  scheme 
that  was  correlated  with  the  surface  observation  data  in  Tables  A6.1  and 
A6.2,  respectively,  for  observation  site  S-1  (Fig.  A 2  and  Table  A2). 
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Elevation  Above  Ground  Level  (m)  Elevotion  Above  Ground  Level  (m) 


Figure  A7.  Project  WIND  stable  (a)  and  unstable  (b)  meteorological  upper  air  initialization 
data  correlated  with  the  surface  observation  data  in  Tables  A7.1  and  A7.2, 
respectively. 


Elevation  Above  Ground  Level  (m)  Elevation  Above  Ground  Level  (m) 


Figure  A8.  Czechoslovakian  stable  (a)  and  unstable  (b)  modeled  upper  air  initialization 
scheme  that  was  correlated  with  the  user  specified  conditions  in  Ttables  A8.1 
and  A8.2,  respectively.  Both  the  stable  and  unstable  profiles  were  applied  to 
observation  sites  S-1  and  S-2  (see  Fig.  A4  and  Table  A4). 
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Table  A12.  Approximation  for  z0  at  the  NTC  from  Arya  (1988), 
and  the  scaling  of  this  z0  value  to  the  normalized  variability  of  the 
Czechoslovakian  site’s  elevation  data  (see  Table  5).  The  Czechoslo¬ 
vakian  site  elevation  data  had  38%  of  the  NTC  site’s  elevation  data 
variability.  The  possible  range  of  z0  for  hilly-mountainous  terrain  was 
set  at  1.93-2.50  (from  Arya,  1988).  The  magnitude  of  this  range  (0.57) 
was  then  scaled  to  the  relative  variability  of  the  Czechoslovakian  site’s 
terrain  elevation  data  (e.g.,  (0.57  x  0.38)  +  1.93  =  2.15). 


Site 

STD. DEV. (elevation) 

Relative  Variability 

raj 

NTC 

172.57 

172.57/172.57  =  1.0 

EH 

Czechoslovakia 

65.76 

65.76/172.57  =  0.38 

2.15 

Table  A13.  The  friction  velocity,  u*,  in  m/s  from  the  10  m  initial¬ 
ization  conditions  for  the  NTC  and  Czechoslovakian  sites  (Tables  A6 
and  A8,  respectively)  using  equation  A-5. 

||  Site  ||  Stable  u*  |j  Unstable  u*  ]{ 


NTC 

Czechoslovakia 


0.2972 

0.2654 


0.8916 

0.7963 


APPENDIX  B 


Sample  of  Analysis  Software  Source  Code 


The  following  source  code,  written  in  a  Kemighan  and  Ritchie  style  C  com¬ 
puter  programming  language,  is  a  sample  of  the  software  that  was  developed  to 
carry  out  the  analysis  of  wind  model  outputs  (discussed  in  the  Results  section). 
Figure  B1  illustrates  the  file  structure  of  this  particular  program  which  was  used 
to  process  100  m  grid  OBWIND  and  WOCSS  results.  In  the  100  m  grid  results, 
each  model  produced  40,401  u  components,  as  well  as  40,401  v  components,  re¬ 
sulting  in  a  total  of  161,604  combined  data  components  per  layer  (e.g.,  at  10,  25, 
and  100  m  above  ground  level). 

The  program  described  here  performs  minimum,  mean,  maximum,  standard 
deviation,  root-mean-squared  difference,  and  agreement  index  calculations  on  the 
separated  u  and  v  components  of  selected  OBWIND  and  WOCSS  wind  vector 
files  for  one  layer  per  program  execution.  The  program  also  solves  for  the  u  and 
v  vector  differences  over  all  grid  points  and  calculates  the  speed  of  any  selected 
u  and  v  component  fields. 

Program  users  specify  the  names  of  the  two  uv  data  files  of  interest,  as  well 
as  the  name  of  the  output  file  near  the  top  of  variables.}*,  which  is  an  include 
file  in  the  main  program,  shell. c.  Each  significant  analysis  step  carried  out  in 
the  program  was  modularized  so  that  users  may  quickly  de-select  unnecessary 
calculations  to  reduce  program  run-time  (which  was  typically  115  seconds  for  a 
complete  analysis  run  on  a  Sun  SPARCstation  IPC  computer).  Output  from  a 
single  execution  of  this  program  would  look  similar  to  the  following  example: 
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WSMRstats. 100m  (Sample  Program  Output  Data  File) 


min  ul  component  ■  -0.462111 
Avg  ul  component  «  -0.575035 
max  ul  component  ■  -0.788063 
min  u2  component  ■  0.000000 
Avg  u2  component  -  -1.190638 
max  u2  component  ■  -1.460600 
ul  std. deviation  ■  0.038507 
u2  std. deviation  ■  0.350084 
rmsd  for  ul-u2  *  0.716975 
ul-u2  agmnt  indx  »  0.999990 


min  vl  component  «  -0.250070 
Avg  vl  component  »  -0.373936 
max  vl  component  ■  -0.502560 
min  v2  component  ■  0.000000 
Avg  v2  component  ■  -0.859763 
max  v2  component  ■  -1.460600 
vl  std. deviation  »  0.050773 
v2  std. deviation  »  0.253497 
rmsd  for  vl-v2  »  0.555462 
vl-v2  agmnt  indx  *  0.999986 
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/a***************#***********************************************/ 

/*  PROGRAM:  SHELL. C;  */ 

/*  Date:  06/05/92;  */ 

/*  Programmer:  Martin  Lee;  *7 

/*  Purpose:  Analyze  and  compare  OBWIND  4  WOCSS  model  results;  •/ 

/*  Use:  sys/file.h  */ 

/*  stdio.h  */ 

/*  math.h;  */ 

/*  variables.h  allows  user  to  specify  file  parameters;  */ 

/*  read  wind  data  into  program  -  read_uv.c;  */ 

/*  a)  uv  wind  file  ■  "wsmr . stable . lvlluv"  (example);  */ 

/*  b)  uv  wind  file  -  "wsmr .unstable. lvlluv" (example) ;  */ 

f*****#*#**************#*****************************************/ 


♦include  <sys/file.h> 

♦include  <stdio.h> 

♦include  <math.h>/*  compile  with  -lm  for  access  to  math  library  */ 
♦include  "variables.h" 

♦include  "diagnostics.c" 

♦include  "read_uv.c" 

♦include  "basic_stats ,c" 

♦include  "rmsdjstats .c" 

♦include  "index. c" 

♦include  "deltas_speed.e" 

♦include  "output. c" 


main  () 

{ 

printf ("\n  MAIN  calling  OBWIND  uv  READ ! \n\n") ; 

read_uv(UVl,ul,  vl) ;  /*  for  OBWIND  u-v  components  */ 

printf ("\n  MAIN  calling  WOCSS  uv  READ ! \n\n") ; 

read_uv ( UV2 , u2 , v2 ) ;  /*  for  WOCSS  u-v  components  */ 

printf ("\n  MAIN  calling  basic  statistics  functions ! \n\n") ; 
basic^statistics <4ulstats,ul) ;  /*  for  OBWIND  u  component  */ 

basic_statistics (svlstats, vl) ;  /*  for  OBWIND  v  component  */ 

basic_3tatistic3 (4u2stats, u2) ;  /*  for  WOCSS  u  component  */ 

basic_statistics(4v2stats,v2) ;  /*  for  WOCSS  v  component  */ 

/*  diagnostic^ <) ;  /*  preliminary  stats  */ 

printf ("\n  MAIN  calling  RMSD  agreement  index  functions! \n\n") ; 
rmsd_stats (ul, vl, u2, v2) ;  /*  root  mean  square  dif.  */ 

printf <"\n  MAIN  calling  index  calculation ! \n\n") ; 
indexO;  /*  agreement  index  */ 

/*  diagnostic^ () ;  /*  final  check  */ 

printf ("\n~MAIN  calling  difference  calculation ! \n\n") ; 
difference!);  /*  u  6  v  subtractions  */ 

printf ("\n  MAIN  calling  speed  calculation!\n\n") ; 

/*  wocss_speed() ;  /*  WOCSS  |uv|  speeds  */ 

printf ("\n  MAIN  calling  output  routine ! \n\n") ; 

output  statsO;  /*  save  analysis  results  */ 

printf“("\n  END  SHELL. C\n") ; 
return (0) ; 


) 

/ 


END  shell. c 


/ 


/**************»******»**«*******• 
/*  variables. h; 

/*  Date:  06/05/92; 

/*  Programmer:  Martir  Lee; 

/*  Purpose:  declare  glcbal  program 

/****••**••*»**,•»*.••„«...*•* 

♦define  X  201 
♦define  Y  201 
♦define  SIZE  40401 
♦define  UV1  "LEVEL3-UV" 

♦define  UV2  "UV3 .DATA" 

♦define  STATS  "WSMRstats . 100m" 


. . *»**•**•«**/ 

*/ 

*/ 

*/ 

variable*.  */ 

. . . 

*  number  of  x  grid  points  */ 

*  number  of  y  grid  points  */ 

*  total  number  of  points!  */ 

/*  OBWIND  100m  UV  DATA  */ 

/*  WOCSS  100m  UV  DATA  */ 

/*  statistical  results  file  / 


♦define  ABS(a)  (f loat) (sqrt < (double) (a) * (a) ) ) 

♦define  SQR(a)  ((a)* (a)) 

♦define  pi  3.14159265 

/A************************** ••**»**********•••*•**••••*••***•* j 


struct  grid_pt  ( 


) ; 


float  component; 


struct  grid_pt  ul [X] [Y] , vl [X] [Y] ; 
struct  gridjpt  u2 [X] [Y] , v2 (X) (YJ ; 


struct  data_atats  ( 
float  avg; 
float  data_avg; 
float  min_data; 
float  max_data; 
double  std; 
double  data_std; 
float  data_threshold; 

) ; 

struct  data_stats  ulstats, vlstats; 
struct  dita_stats  u2stats, v2stat.s; 

/****************!»*************#***********#******************/ 
FILE  *fp,  *u_deltas, *v_deltas, *uv_speed, "output; 

. . . . . . . . . 


char  filename[81] ; 
int  i,j,num; 

float  value, u_d iff, v_di f f , speed ; 

double  u_register, v_register, delta_u,delta_v, rmad_u, rmsd_v; 
double  inner_u_sura,  w_u__prime,  o_u__prime,  inner_v_sum; 
double  w_v_p r ixne ,  o_v_p r  ime ,  d_u ,  d_v  ; 

/********************  END  variables. h  ***********************/ 
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/*  diagnostics .c  •/ 
/*  Date:  06/05/92;  ./ 
/*  Programmer:  Martin  Lea;  »/ 
/*  Purpose:  error  detection/trapping;  »/ 
/* . ***************************** . * . * . —  / 


void 

diagnostic_l {) 

{  /*  BEGIN  first  diagnostic  */ 

printf ("\n  in  first  diagnostic ! \n\n") ; 
printf ("Average  ul  component  -  %f\n",ulstats.data_avg) ; 
printf ("Average  vl  component  -  %f \n", vlstats .data_avg) ; 
printf ("min  ul  component  -  %f\n", ulstats .min_data) ; 
printf ("max  ul  component  *  %f \n", ulstats .max^data) ; 
printf ("min  vl  component  *  %f\n", vlstats .min_data) ; 
printf ("max  vl  component  *  %f\n\n", vlstats ,max_data) ; 
printf ("Average  u2  component  -  %f \n", u2stats .data_avg) ; 
printf ("Average  v2  component  -  %f \n", v2stats .data_avg) ; 
printf ("min  u2  component  -  %f \n",  u2stats .min_data) ; 
printf ("max  u2  component  -  %£\n",  v2stats.max_data) ; 
printf ("min  v2  component  -  %f\n",u2stats.min__data) ; 
printf  ("max  v2  component  •  %f \n",  v2stats  .max__data) ; 
printf  ("\n  exiting  first  diagnostic ! \n\n") ; 

}  /*  END  first  diagnostic  */ 


void 

diagnostic_2 ( ) 

{  /*  BEGIN  second  diagnostic  */ 

printf ("\n  in  second  diagnostic ! \n\n") ; 
printf ("i-  %d,  j-  %d\n",i,j); 
printf ("ul-  1f\n",ul (i] [ j] .component) ; 
printf ("vl-  %f \n", vl [ij ( j] . component) ; 
printf ("u2-  %f\n",u2 [i] ( j] .component) ; 
printf ("v2-  %f\n", v2 [i] t j] .component) ; 
printf ("u_register  -  %f \n", u_register) ; 
printf ("v_register  -  %f \n"f v_register) ; 
printf  ("\n  exiting  second  diagnostic ! \n\n") ; 

)  /*  END  second  diagnostic  */ 


void 

diagnostic_3 {) 

(  /*  BEGIN  third  diagnostic  */ 

printf ("\n  in  third  diagnostic ! \n\n") ; 
printf ("rmsd  for  u  -  %f\n", (float) rmsd_u) ; 
printf ("rmsd  for  v  -  %f\n", (float) rmsd_v) ; 
printf ("agreement  index  for  u  -  %f\n",d_u); 
printf ("agreement  index  for  v  -  %f\n",d_v); 
printf  ("\n  exiting  third  diagnostic ! \n\n") ; 

)  /*  END  third  diagnostic  */ 

END  diagnostics .c  ********** 


HI 


/ 


/ 


. . . . . ****** . ***•*/ 

/*  read  uv.c;  */ 
/*  Date?  06/05/92;  */ 
/*  Programmer:  Martin  Lee;  */ 
/*  Purpose:  Read  UV  wind  data  files;  */ 
. . * . * . * . * . **•**/ 


void 

read_uv (uv_f ile, u, v) 
char  uv_file{]; 
struct  gridjpt  ulX)(Y); 
struct  gridjpt  v[X)(Y]; 
< 


strcpy (filename, uv_f ile) ; 
fp  -  fopen {filename,  "r") ; 

if  (fp  —  NULL) 

{ 

printf ("Error  on  open\n"); 
exit  ( 0 ) ; 

) 

rewind (fp) ; 

for  (i-0;  i<«X-l;  i++) 

I 

for  (j»0;  j<-Y-l;  j++) 

( 

num  -  fscanf (fp, "%f", Svalue) ; 
if (num  !-  1) 

( 

printf ("Error  on  u  read\n"); 
exit (0) ; 

) 

u(i] tj] .component  -  value; 

} 


for  (i-0;  i<-X-l;  i++) 

{ 

for  (j-0;  j<-Y-l;  j++) 

< 

num  -  fscanf (fp, "If ", (value) ; 
if (num  !-  1) 

( 

printf ("Error  on  v  read\n") ; 
exit ( 0 ) ; 

) 

v(i] ( j) .component-  value; 

) 


) 

/*****«*«****••****•••**  end  read  uv.c  ************************/ 
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Z**************************************************. »»*.*. *••**.* 

/*  basic  stats. c 
/*  DATE  :~0 6/ 05/92 ; 

/*  PURPOSE:  perform  basic  statistical  tests  on  input  data  array; 
/*  PROGRAMMER:  Martin  Lee; 

/*•»*»*******»****•******»•»•***•*•*»*****•«•*•*****. ............ 

void 

basic_statistics (stats, data) 
struct  data_stats  *stats; 
struct  grid_pt  data{X](Y]; 
l 


*/ 

*/ 

•/ 

*/ 

*/ 

*/ 


/*  BEGIN  Calculation  of  the  MEAN  data[X][Y]  components  */ 
stat3->avg  -  0.0; 
for  (i-0;  i<X-l;  i++> 

{ 

for  (j-0;  j<Y-l;  j++> 

stats->avg  -  stats->avg  +  data (i] ( j] .component; 

) 

stats->data_avg  -  stats->avg/SZZE; 

/*  END  Calculation  of  the  MEAN  d(X] [Y]  components  */ 

/*  BEGIN  Calculation  of  the  MIN,  and  MAX  data[X] [Y]  components  */ 
3tats->min_data-3tats->data__avg; 
stats->max  data-stats->data  avg; 
for  (i-0;  I<X-1;  i++) 

( 

for  (j-0;  j<Y-l;  j++) 

( 

if (ABS (data[i] [j] .component)  <  A8S (stats->min_data) ) 
3tats->min_data  -  data [ i ][ 1 1 . component ; 
if (ABS (data [i] (j] .component)  -  ABS (stats->max_data) ) 
stats->max  data  -  datafi] [ j] .component; 

) 

) 

/*  END  Calculation  of  the  MIN,  and  MAX  data[X][Y)  components  */ 

/*  BEGIN  Standard  Deviations  of  the  data[X][Y]  components  */ 

3tats->std  -  0.0; 
for  (i-0;  KX-l;  !+■*•) 

( 

for  (j-0;  j<Y-l;  j++) 

3tats->std  -  stats->std+Ss3R(data(i] l jj .component-stat3->data  avg); 

) 

atats->data_3td  -  sqrt (stats->std/SIZE) ; 

/*  END  Standard  Deviations  of  the  data[X][Y]  components  */ 

/*  BEGIN  find  | thresholds |  of  the  data[X][Y]  components  */ 

stats->data_threshold  -  ABS (stats->dsta_avg)  +  stats->data__std; 

/*  END  find  upper  thresholds  of  the  data[X][Y]  components  */ 


J 

/*«•***•*•**•*•***«•*  end  basic  stats. c  ****t*»*o*****«**ot»/ 
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/*  rmsd_stats .  c  */ 
/*  DATE:  06/05/92;  */ 
/*  PURPOSE:  carry  out  RMSD  calculation  on  wind  modal  data;  •/ 
/*  PROGRAMMER:  Martin  Lea;  */ 
. . . . . ***** . . . 


void 

rmad_stat3 (ulr, vlr, u2r, v2r) 
struct  gridjpt  ulr[X][Y]; 
struct  grid_pt  vlr[X][YJ; 
struct  grid_pt  u2r(X][Y]; 
struct  grid_pt  v2r{X][Y]; 

{ 

/*  BEGIN  RMSD  Calculation  for  ulr-vlr,  and  u2r-v2r  wind  fields  */ 

delta_u  -  0.0; 
delta_v  -  0.0; 
delta_u  -  0.0; 
delta_v  -  0.0; 

for  (i-0;  i<-X-l;  i++) 

( 

for  (j-0;  j<-Y-l;  j++) 

1 


delta_u  -  (u2r(i] tj] .component)- (ulr [il f j ) .component) ; 
u_register  -  SQR(delta_u)  +  u_register; 


delta_v  -  (v2r[i] £  j  3 .component) - (vlrli] [j] .component) ; 
v_register  ■  SQR(delta_v)  +  v_register; 

/*  diagnostic _2();  /*  rmsd  error  trap  */ 

) 

) 

rmsdju  ■  sqrt (u_register/SIZE) ; 
rroad__v  -  sqrt  (v_register/SIZE) ; 

/*  END  RMSD  Calculation  for  ulr-vlr,  and  u2r-v2r  wind  fields  */ 

I 

/********************  end  rmsd  stats. c  **********************/ 
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. . . . . . 

/*  index. c  •/ 

i*  DATE:  06/05/92;  */ 

/*  PURPOSE:  agreement  index  test  of  wind  model  data;  •/ 

/*  PROGRAMMER:  Martin  Lee;  */ 

/.ft**.*******.**************...... .......ft**.*******. 


void 
index ( ) 

( 

/•  BEGIN  agreement  index  of  ulr-vlr,  and  u2r-v2r  wind  fields  */ 

inner_u__sum  -  0.0; 
w_u_prime  •  0.0; 
o_u_prime  »  0.0; 
inner_v_sum  “  0.0; 

w_v_p r ime  *  0.0; 
o_vjprime  -  0.0; 

for  (i-0;  i<“X-l;  i++) 

{ 

for  (j-0;  j<«Y-l;  j++) 
i 

w_u_prime  -  ABS (u2 (i] t j] .component  -  ulstats .data_avg) 
o__u_jprime  -  ABS (ul [i]  { j]  .component  -  ulstats .data_avg) 
inner_u_aum  -  inner_u_3um  +•  SQR(w__u_prime  +  o_ujprime) 


w__v_jprime  »  ABS <v2 (ij  (j]  .component  -  vlstat3.data_avg) 
o__vjprime  -  ABS (vl [ij (j] .component  -  vlatat3.data_avg) 
inner_v_sum  -  inner_v_aura  +  SQR(w_v_prime  +  o_vjprime) 

) 

} 

d_u  -  1.0  -  SIZE*SQR(rm3d_u) /inner_u_3um; 
d_v  -  1.0  -  SIZE*SQR(rmsd_v) /inner_v_sura; 

/*  END  agreement  index  of  ulr-vlr,  and  u2r-v2r  wind  fields  */ 

> 

/*•■»•#*******»******•**  end  index. c  ************************/ 
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/*  deltas_speed.c 
/*  DATE:  06/05/92; 

/*  PURPOSE:  find  differences  and  speeds  of  u  and  v  fields; 
/*  PROGRAMMER:  Martin  Lee; 


void 

difference {) 

( 

if  ( (u_deltas-f open  ("u. deltas”,  "w+") )  —NULL) 


< 


fprintf (stderr, "cannot  open  u. deltas !\n") ; 
exit (1) ; 

) 

rewind (u_delt as) ; 

if  ( (v_deltas-f open  ("v. deltas",  "w+") )  —NULL) 

I 

fprintf (stderr, "cannot  open  v.deltas!\n") ; 
exit (1) ; 

) 

rewind (v_delt as) ; 
for  (i-0;  i<-X-l;  i++) 

{ 


for  (j-0;  j<-Y-l;  j++) 


{ 


) 


'*/ 

*/ 

*/ 

*/ 

•/ 

,*/ 


u_diff  -  u2 [i] [j] .component  -  ul [i] { j] .component; 
fprintf (u_deltas, "%f \n",u_diff ) ; 

v_diff  -  v2 [i] [ j] .component  -  vl[i] [j] .component; 
fprintf (v_deltas, "%f \n", v_diff ) ; 


/*  end  differences  */ 


void 

wocss_speed() 

{ 

if  ( (uv  speed- f open  ("uv. speed",  "w+") )  —NULL) 

< 

fprintf (stderr, "cannot  open  u. deltas ! \n" ) ; 
exit (1) ; 

) 

rewind (uv_s  peed) ; 
for  (i-0;  K-X-l;  i++) 

( 

for  (j-0;  j<-Y-l;  j++) 

( 

speed  -  sqrt (SQR(u2 [i] [ j] .component) +SQR (v2 fij [j] .component) ) ; 
fprintf (uv_speed, "%f \n", speed) ; 

) 

) 

}  /*  end  wocss_speed  */ 

/********************  end  deltas_speed.c  ***********************/ 
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/ *  output . c  *  / 
/*  DATS:  06/05/92;  */ 
/*  PURPOSE:  file  analysis  results;  */ 
/*  PROGRAMMER:  Martin  Lee;  V 


/*******#************************★******************«*******«*/ 

void 

output  stats <) 

{ 

st rcpy (filename, STATS) ; 
output  ■  f open (filename,  "w+"); 

if  (output  »■  NULL) 

( 

printf ("Error  on  open\n") ; 
exit (0) ; 

) 

rewind (output) ; 

fprintf (output, "%s\n", filename) ; 
fprintf (output, "\n") ; 

fprintf (output, "min  ul  component  -  %f\n", ulstats .min_data) ; 
fprintf (output, "min  vl  component  -  %f\n", vlstats .min_data) ; 
fprintf (output, "Avg  ul  component  -  %f\n",  ulstats. data_avg) ; 
fprintf (output, "Avg  vl  component  *  %f \n",  vlstats .data_avg) ; 
fprintf (output, "max  ul  component  -  %f \n", ulstats .max_data) ; 
fprintf (output, "max  vl  component  ■  % f\n", vlstats. max_data ) ; 
fprintf (output, "\n") ; 

fprintf (output, "min  u2  component  ■  %f \n", u2stats .min_data) ; 
fprintf (output, "min  v2  component  •  %f\n*,u2stats.min_data) ; 
fprintf (output, "Avg  u2  component  -  %f \n", u2stats ,data_avg) ; 
fprintf (output, "Avg  v2  component  -  %f \n", v2stats ,data_avg) ; 
fprintf (output, "max  u2  conponent  -  *f \n", v2stats .max_data) ; 
fprintf (output, "max  v2  component  ■  %f\n", v2stats .max_data) ; 
fprintf (output, "\n") ; 

fprintf (output, "ul  std. deviation  -  %f\n", ulstats .data_std) ; 
fprintf (output, "vl  std. deviation  -  %f\n", vlstats. data_std) ; 
fprintf (output, "u2  std. deviation  -  %f \n", u2stats ,data_std) ; 
fprintf (output, "v2  std. deviation  •  %f\n", v2stats.data_std> ; 
fprintf (output, "\n") ; 

fprintf (output, "rmsd  for  u  -  %f\n", (float) rmsd_u) ; 
fprintf (output, "rmsd  for  v  -  lf\n", (float) rmsd_v) ; 
fprintf (output, "\n") ; 

fprintf (output, "agreement  index  for  u  -  %f\n",d_u); 
fprintf (output, "agreement  index  for  v  ■  %f\n",d_v); 


} 

/ 


END  output. C 


/ 


APPENDIX  C 

Summary  of  Software  Execution  Details 

Both  the  OBWIND  model  (841  lines  of  source  code)  and  the  WOCSS  model 
(4,343  lines  of  source  code)  were  installed,  configured,  compiled,  and  tested  prior 
to  commencement  of  the  model  runs.  Each  model  was  configured  to  produce 
windfield  outputs  at  10,  25,  and  100  m  above  ground  level.  Each  model  was 
then  sequentially  initialized  and  executed  with  the  correct  terrain  and  meteoro¬ 
logical  data  for  each  study  site.  Figure  Cl  maps  out  all  of  the  OBWIND  and 
WOCSS  model  runs  that  were  executed  to  satisfy  the  cross-model  study  design 
goals  (illustrated  in  Fig.  10).  Execution  of  all  OBWIND  and  WOCSS  model 
runs  indicated  in  Fig.  Cl  produced  5,688,360  wind  field  data  elements.  These 
wind  field  data  elements  represented  the  orthogonal  x  and  y  axes  components 
of  individual  wind  vectors,  V,  produced  by  the  model  runs  at  each  grid  point 
(e.g.,  V  =  tii  +  vj,  where  u  =  ui  and  v  =  v j  are  the  x  and  y  axes  component 
vectors,  respectively;  and  i  and  j  are  the  respective  unit  vectors  along  the  x  and 
y  axes).  Table  Cl  summarizes  the  number  of  u  and  v  wind  field  vector  compo¬ 
nents  produced  during  each  model  run,  and  Table  C2  summarizes  the  cumulative 
total  of  these  components  produced  from  all  the  model  runs  executed.  Sizes  of 
the  resulting  computer  files  containing  these  wind  field  components  are  noted  in 
Table  C3. 

Table  Cl.  Number  of  individual  wind  component  data  elements  pro¬ 
duced  during  both  an  OBWIND  and  WOCSS  model  run. 


Grid  Spacing  (m) 

u  components 

v  components 

1000  (each  model  run/per  layer) 

441 

441 

250  (each  model  run/per  layer) 

6,561 

6,561 

100  (each  model  run /per  layer) 

40,401 

40,401 

Total  output/model/layer  [  47,403 

47,403 

|  Combined  outputs/model/layer  | 

04,1 

906 
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Model  Run  Caws  Executed  using  both  OOWIND  and  WOCSS 


Figure  Cl .  OBWIND  and  WOCSS  model  runs  were  completed,  and  the  horizontal  uv  wind  vector  data  files 
were  saved  lor  10, 25,  and  100  m  above  ground  level  in  the  file  directory  scheme  illustrated  above. 


Table  C2.  Total  number  of  individual  data  elements  produced  during 
OBWIND  and  WOCSS  model  runs  for  WSMR,  NTC,  Project  WIND, 


and  Czechoslovakia. 


Model  Runs 

u  components 

v  components 

OBWIND  +  WOCSS  output/layer 

94,806 

94,806 

Stable  -f  Unstable  Cases/layer 

189,612 

189,612 

(WSMR+NTC+Proj.WIND)/layer 

568,836 

568,836 

Czech,  sta.  1  stable/layer 

189,612 

189,612 

Czech,  sta.  2  unstable/layer 

189,612 

189,612 

Total  for  Czechoslovakian  runs 

379,224 

379,224 

Cumulative  Total  for  1  level 

948,060 

948,060 

Cumulative  Total  for  3  levels 

2,844,180 

2,844,180 

u  &  v  elements  from  all  runs 

5,688,360 

Table  C3.  Combined  u  and  v  windfield  output  data  file  sizes  (in 


kilobytes)  for  both  OBWIND  and  WOCSS  model  runs. 


Model  Runs 

100m  grid 

250m  grid 

1000m  grid 

OBWIND/layer/run 

329 

139 

9.3 

WOCSS/layer/run 

329 

139 

9.3 

OBWIND+WOCSS/layer 

658 

278 

18.6 

OBWIND+WOCSS/3  layers 

1,974 

834 

55.8 

Stable(S)+Unstable(U)  runs 

3,948.0 

1,668.0 

111.6 

Cumulative  total  for  S+U  runs 

5,727.6 

Cumulative  toted  for  all  runs 

34,365.6 

To  produce  this  resultant  output  data,  model  runs  required  the  equivalent  of 
4.3  hours  of  uninterrupted,  dedicated  CPU  time  (detailed  in  Tables  C4  and  C5)  on 
a  Sun  Microsystems  SPARCstation  IPC  computer  workstation.  1  The  execution 
times  presented  in  Tables  C4  and  C5  are  very  dependent  on  the  number  of  layers 
solved  for  by  each  model,  as  well  as  on  all  other  user  specified  input/output 
configuration  specifications. 

lA  17.4  million  instructions  per  second  (MIPS),  2.1  million  floating  point  operations  per 
second  (MFLOPS),  25-MHs  SPARC  proccessor  platform  with  8  MegaBytes  of  random  access 
memory  (RAM)  and  a  1.2  GigaByte  hard  disk. 
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Tables  C4  and  C5  imply  that  the  WOCSS  model  executes,  on  an  average 
basis  (over  100,  250,  and  1000  m  grids),  approximately  3  times  faster  than  the 
OBWIND  model.  The  reduction  in  WOCSS  execution  time  is  primarily  due  to 
differences  between  the  OBWIND  and  WOCSS  models’  interpolation  schemes 
and  internal  algorithm  process  implementation. 


Table  C4.  Typical,  dedicated  CPU  time  (in  seconds)  required  to  make 
OBWIND  and  WOCSS  model  runs  -  where  all  standard  output  was 
redirected  to  a  file  (e.g.,  no  screen  outputs  were  presented  during 
model  execution  to  reduce  bias  in  program  run-times). 


Grid  Spacing  (m) 

OBWIND 

WOCSS 

100  (per  run) 

994 

327 

250  (per  run) 

155 

61 

5.5 

5.4 

Total  (per  model  run) 

1,154.5 

393.4 

Combined  Total/run 

1,547.9 

Table  C5.  Total  estimate  of  dedicated  CPU  time  (in  seconds)  required 
to  make  all  cross-case  OBWIND  and  WOCSS  model  runs  -  where 
standard  output  was  redirected  to  a  file  (e.g.,  no  screen  outputs  were 
presented  during  model  execution  to  reduce  bias  in  program  run¬ 
times). 


|  Model  runs 

OBWIND 

WOCSS  1 

Total  WSMR  runs 

2,309.0 

786.8 

Total  NTC  runs 

2,309.0 

786.8 

Total  Proj.WIND  runs 

2,309.0 

786.8 

Total  Czech,  runs 

4,618.0 

1,573.6 

Cumulative  Total 

11,545.0 

3,934.0  | 

Combined  Cum.  Total 

15,479.0 
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